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High temperature superconductivity in iron pnictides and chalcogenides emerges when a magnetic 
phase is suppressed. The multi-orbital character and the strength of correlations underlie this 
complex phenomenology, involving magnetic softness and anisotropies, with Hand’s coupling playing 
an important role. We review here the different theoretical approaches used to describe the magnetic 
interactions in these systems. We show that taking into account the orbital degree of freedom allows 
us to unify in a single phase diagram the main mechanisms proposed to explain the (tt, 0) order in 
iron pnictides: the nesting-driven, the exchange between localized spins, and the Hand induced 
magnetic state with orbital differentiation. Comparison of theoretical estimates and experimental 
results helps locate the Fe superconductors in the phase diagram. In addition, orbital physics is 
crucial to address the magnetic softness, the doping dependent properties, and the anisotropies. 


I. INTRODUCTION 

The discovery of iron superconductors mm challenged 
the uniqueness of cuprates as high Tc superconductors, 
triggering the search for new materials in previously un¬ 
explored directions. Most iron superconductors are mag¬ 
netic when undoped and the suppression of the mag¬ 
netism via doping, pressure or isoelectronic substitution 
is accompanied by the appearance of a superconducting 
phase I3H5]. There is also a structural phase transition 
that may coincide with the magnetic transition or occur 
at a higher temperature. Details may depend on the par¬ 
ticular structure or chemical composition but the general 
trends are robust. All the families of iron superconduc¬ 
tors have atomic layers with Fe in a square lattice and a 
pnictogen (P, As) or a chalcogen (Se, Te) out of plane in 
tetrahedral coordination [6], see Fig. Understanding 
the magnetism and electronic correlations in these sys¬ 
tems may be the clue for elucidating the yet unknown 
pairing mechanism for high Tc superconductivity mm- 

Superficially, the iron superconductors phase diagram 
resembles very much that of the cuprates in the sense that 
a magnetic phase is replaced by a superconducting one 
with doping [8]. However, there are some deep differences 
whose relevance for superconductivity is still to be deter¬ 
mined. Undoped cuprates are well known for being Mott 
insulators and (7r,7r) order antiferromagnets. On the 
other hand, undoped Fe superconductors are generally 
metallic and show different types of magnetic arrange¬ 
ments The one that dominates is the columnar 

(tt, 0) one (with axis defined in the Fe-Fe first neighbours 
direction): antiferromagnetic in the x-direction and fer¬ 
romagnetic in the ^-direction, see Fig. This order (or 
related magnetic fluctuations) is usually accompanied by 
a tetragonal to orthorhombic structural transition and 
in-plane {x/y) anisotropies. 
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The (tt, 0) order is found in the undoped iron arsenides, 
a denomination that includes different chemical composi¬ 
tions and structures -all with FeAs planes. The arsenides 
mainly comprise the 1111 family, like LaOFeAs, the 122, 
like BaFe 2 As 2 , and the 111, like LiFeAs and NaFeAs. 
The related phosphides (with FeP planes) are non mag¬ 
netic and their superconducting TcS are low [T]. Most 
of the arsenides show the columnar order with a small 
magnetic moment, with typical values under I^b per Fe, 
and are metals [9] . LiFeAs is an exception to the rule as 
it does not order magnetically, it does not go through a 
structural transition and it is superconductor [10]. The 
iron chalcogenide FeTe has a different magnetic order, 
the double stripe, namely a double FM column along the 
diagonals (the Fe-Fe second neighbours direction) and a 
larger magnetic moment ^ 2ys per Fe m- On the other 
hand, bulk FeSe suffers a structural transition at a much 
higher temperature than the superconducting one but 
blue static magnetism only arises under pressure nansj. 
Related systems arise when some spacer is introduced be¬ 
tween the FeSe layers [min]. The parent compounds of 
the alkaline Fe selenides A^Fe 2 -a:Se 2 (A=K, Rb, Cs, Tl), 
with Fe vacancies arranged in a particular pattern, have 
a block antiferromagnetic order with a large magnetic 
moment of ^ S/i^ per Fe, and are insulators Enmn. 
Transitions between different magnetic orders have been 
observed as a function of doping [22] and pressure [23ll^ . 

The electronic structure gives us hints to under¬ 
stand the observed metallicity and the differences with 
cuprates. Models for cuprates usually involve a single 
orbital. The situation is very different in the iron super¬ 
conductors. Density Functional Theory (DFT) calcula¬ 
tions show that the density of states around the Fermi 
level is dominated by the Fe-pnictogen (or chalcogen) 
planes [251127] . In undoped compounds, Fe is in a 3d^ 
valence state and the crystal field is much smaller than 
the bandwidth. This implies that all the d orbitals may 
play a role at low energies. The Fermi surface of iron 
superconductors has both electron and hole pockets at 
different /c-points and different orbitals dominate in dif¬ 
ferent parts of the Brillouin zone, see Fig. There¬ 
fore, multi-orbital models are required to describe the Fe 
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FIG. 1: (Left) The crystal structure of all the families of Fe 
superconductors share planes with Fe in a square lattice with 
a pnictogen or a chalcogen (here As) located out of plane in 
a tetrahedral coordination. The relative position of the out- 
of-plane anions with respect to the plane can be described 
by the angle a. Reproduced from Ref. [ST]. (Right) Typical 
Fermi surface defined in the unfolded Brillouin zone (which 
corresponds to a single Fe unit cell with the x and y direc¬ 
tions defined along the Fe nearest neighbours). This Fermi 
surface has been calculated with the tight-binding proposed 
in Ref. [31] with the a corresponding to a regular tetrahedron. 
The arrow indicates the wave-vector of the columnar ordering 
(tt, 0). The pockets at F and M are hole pockets and the ones 
at X and Y are electron pockets. The color code refers to 
the orbital with the highest weight on a particular part of the 
Fermi surface: Blue for xy, red for zx and green for yz. The 
eg orbitals 3z^ — and — y‘^ are also present at the Fermi 
surface though with a smaller weight [ST]. 
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superconductors. Although most of these materials are 
metallic when undoped, correlations are known to play 
an important role as, comparing ab-initio bands with An¬ 
gle Resolved Photo Emission Spectroscopy (ARPES) and 
quantum oscillations measurements, a band renormalisa¬ 
tion of ^ 3 has been estimated [28H3Q] . 

The first attempts to understand the phase diagrams 
of the Ee superconductors came from two opposite ap¬ 
proaches: the weak and the strong coupling limits. The 
weak coupling approach was supported by the metallic- 
ity of the parent compounds, the low magnetic moment, 
and the apparent nesting on the Eermi surface with a 
Q = (7i',o) wave-vector EZllSl], see Eig. The strong 
coupling approach |33l [34] was motivated by the observed 
band renormalisation, the bad metallicity and later by 
the discovery of families of Ee superconductors whose 
magnetic order could not be explained by Eermi surface 
nesting. The multi-orbital character of the Ee supercon¬ 
ductors has inspired models in which itinerant and lo¬ 
calised electrons coexist [35H38] . similarly to other tran¬ 
sition metal oxides. Theoretical studies of multi-orbital 
systems found the possibility of having an orbital selec¬ 
tive Mott transition [391441] . But the complexity did not 
end here, and the effect of correlations in multi-orbital 
systems has become a research field on its own right. 
Hund’s coupling is one of the main characters of the emer¬ 
gent story, leading to a correlated Hund’s metal, jiaiis] 
with a suppressed but non-zero quasiparticle weight. Dif¬ 
ferent degrees of correlation for the orbitals (orbital dif¬ 
ferentiation), as observed in ARPES experiments [41], 
adds up to this richness. 


An additional complexity arises in connection with the 
anisotropies observed in these materials. The (tt, 0) mag¬ 
netic order and the tetragonal to orthorhombic structural 
transitions break the four-fold rotational symmetry of the 
lattice. The two transitions usually happen very close in 
temperature revealing a strong connection between the 
spin and lattice degrees of freedom. The in-plane resis¬ 
tivity is anisotropic in an unexpected way as it is larger 
for the shorter ferromagnetic direction. In addition, the 
magnitude of the anisotropy is largest when the struc¬ 
tural distortion and the magnetism are weaker [444148] . 
Signatures of anisotropy have also been found in other 
experiments [494(65] . The degeneracy of yz and zx is 
broken in both the magnetic and non magnetically or¬ 
dered phases below the structural transition as revealed 
by ARPES [664(68] and X-ray experiments [69]. A ne¬ 
matic phase, characterised by the reduction of the point 
group symmetry from C 4 tetragonal to C 2 orthorhom¬ 
bic, is found between the magnetic and structural tran¬ 
sitions. Nematic transitions are also believed to occur in 
Sr 3 Ru 2 07 and cuprates CO]. 

There is nowadays an intense debate on the origin of 
the nematicity in the Ee superconductors. Different ex¬ 
periments seem to indicate that it is electronic in ori¬ 
gin [4541471 [ 5 TI |54| [55l |67] but it is difficult to pinpoint 
between the spin miizii or orbital |Ml[I3IZlj degrees 
of freedom due to the spin-orbital entanglement. A re¬ 
cent effective model sensitive to the orbital degree of free¬ 
dom finds that spin fluctuations generate orbital split- 
ting [75|. The possible role of impurities is also un¬ 
der investigation @3 E3 E3 CSHZni. Importantly, or¬ 
bital or spin fluctuations are also likely candidates for 
the superconducting pairing glue: magnetic fluctuations 
are believed to mediate -superconductivity [271180], 
and orbital fluctuations are thought to lead to 5 ++- 
superconductivity [zaEi]. 

In this paper we review the magnetic interactions in Ee 
superconductors with main emphasis on the connection 
between the orbital degree of freedom and the correla¬ 
tions, anisotropies and magnetic softness. In Sec. [n|the 
different models and techniques used for studying the 
magnetic properties of the Ee superconductors are intro¬ 
duced. The following section deals with correlations and 
the (tt, 0) order and is divid ed in f our subsections: weak 
coupling descriptions (Sec. Ill A), strong coupling de¬ 
scriptions (Sec. IIIB| ), ab-initio insights (Sec. IIIC), and 
the role of Hund ’s coupling and the orbital degree of free¬ 
dom (Sec. Sec. [previews the origin of anisotropic 

properties of Ee-superconductors. Sec. |V| deals with the 
current understanding of the observed magnetic softness 
in these systems. We end in Sec. jVT] with the conclusions. 


II. MODELS AND TECHNIQUES 

Magnetic interactions in iron superconductors have 
been studied using multiple techniques. Among them, 
ab-initio methods have played an important role in set- 
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FIG. 2: Some of the different antiferromagnetic orders within 
the Fe plane that appear in Fe-based superconductors. Each 
arrow represents the magnetic moment on an Fe ion. (a) The 
Q = (tt, tt) order, (b) the columnar order Q = (tt, 0), the most 
common order in the Fe-pnictides, and (c) the double stripe 
order as in FeTe. 

ting the essential features of the band structure and 
identifying the magnetic ground state [26l [35l I82H88] . 
DFT methods are single-particle mean field descriptions 
and deviations from their predictions are expected when 
many-body effects play an important role in the non¬ 
magnetic state. 

A different approach deals with lattice multi-orbital 
systems adding interactions to tight-binding models 
which mimic the bandstructure close to the Fermi sur¬ 
face iniEziiHsiiHmi]- Models from 2 to 8 orbitals per 
one-Fe unit cell have been used. Pnictogen and chalco- 
gen orbitals lie several eV below the Fermi level [93] and 
many models only include their effect in the hopping inte¬ 
grals between Fe-d orbitals. Some authors have used only 
2 to 4 orbitals per Fe atom. However the 5 Fe-orbitals 
{yz^ zx, xy, 3z^ —r^ and x^ —y^) are required to account 
in the same model for the Fermi surface, the hoppings 
and the orbital fillings, all of which could be relevant in 
setting the magnetic ground state and correlations. In 
the following we focus in calculations including at least 
the 5 d-orbitals. We use a one-Fe unit cell with x and y 
axis along the nearest Fe-Fe bonds [271|3lll9H96]. 

Typical hamiltonians include the tight-binding, calcu¬ 
lated within Slater-Koster m or by fitting to an ab-initio 
band structure [ssiiHaiiT], the crystal field splitting and 
local interactions restricted to Fe orbitals: intraorbital 
[/, interorbital t/'. Hand’s coupling Jh and pair-hopping 
J'. 

^ 4,'y,crCj,l3,a + h.C. + U 

hi,7^/3,CT i,7 

+ (U' - ^) E - 2Jh E 

i,7//3 i,7,cr 

i,jf label the Fe sites, a the spin and 7,^5 the orbitals. 
Most of the models assume orbital-independent interac¬ 
tions and the relations U' = U — 27//, and J' = Jh m, 
valid in rotationally invariant systems, leaving only two 
independent parameters, U and Jh- Repulsion between 
electrons requires Jh < U/3. Interactions have been esti¬ 
mated with constrained RPA [98], constrained LDA [39| 


or GW [99] methods. For five-orbital models they give 
U ^ 2.5 — 3 eV and Jh ^ 0.35 — 0.45 eV in arsenides and 
U ^ 3 — 4 eV and Jh ^ 0.45 — 0.5 in chalcogenides [98] . 
Models which include p orbitals, as those used in some 
LDA+DMFT [lOOi llOlj or Gutzwiller [102] approaches, 
generally assume non-interacting p orbitals (35] [98] and 
use larger values of U ^4 — 5 eV and Jh ^ 0.5 — 0.8 for 
d-orbitals, as screening from p-orbitals is not included. 

Techniques used to study these multi-orbital mod¬ 
els include: Random Phase Approximation (RPA) [83l 
[T03] . FLEX [91], momentum and real-space Hartree- 
Fock (HF) [37[ |l2] ll04H108j . Functional Renormalisa¬ 
tion Group (FRG) [109] . Dynamical Mean Field The¬ 
ory (DMFT) [^[4^ [T00l[TT0m^ . Monte Carlo ttni, 
Gutzwiller variational theory [TMl rriMiTTl and Slave 
Spin uni nm iTTfl] . Different techniques approach the 
model in weakly or strongly correlated limits and can 
underestimate or overestimate the effect of interactions. 
For this reason and to better understand the underlying 
physics, the phase diagram as a function of U and Jh is 
often explored. 

More simplified descriptions |12QI4125| consider band 
models with electron pockets dX X = (tt, 0) and Y = 
(0,7r) and one or two hole pockets at F = (0,0). Bands 
are frequently assumed to be parabolic with no orbital 
content. Among all possible fermionic interactions only 
density-density and pair-hopping between electron and 
hole pockets, related to the spin-density wave instability 
are kept. These models have been treated with Renor¬ 
malisation Group (RG) [1201112111125] or within mean 
field approaches [12211123] for the collective magnetic de¬ 
grees of freedom. These low energy models are suitable to 
deal with itinerant electrons in weakly interacting met¬ 
als, i.e. when the ordered states are well described by an 
instability of the Fermi surface and if the orbital content 
of the Fermi surface does not play a relevant role. Very 
recently a low energy model, which only keeps interac¬ 
tions in the spin-channel but retains the orbital character 
of the bands and the information about U and Hand’s 
coupling, has been derived Ca¬ 
in the opposite limit, strong coupling approaches [331 
13311371 rr^ ITTt] assume that the electrons are localised 
and describe the system with interacting lattice spin 
models as the Heisenberg model 

= XY ^ E (2) 

' ' (hi) ' ' ((hi)) 

Here Si = is the atomic moment and (i, j) and 

{{i^j)) are restricted to first and second nearest neigh¬ 
bours respectively. The small crystal field splitting in 
iron superconductors suggests |5'| = 2 due to Hand’s rule. 
Ji and J 2 have been estimated applying perturbation 
theory to multi-orbital models [sain and in ab-initio 
calculations [84] [85] 11281 I129j . but in most cases they 
are assumed to be unknown. In a tetragonal state, Ji is 
equal in x and y directions. To fit some neutron mea¬ 
surements, the spin interactions to first nearest neigh- 
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hours are allowed to be anisotropic, i.e. Ji^a 7 ^ 

Some measurements and magnetic orderings have been 
described including in the model a J 3 interaction to third 
neighbours [90] or a biquadratic term K{SiSj)^ to near¬ 
est neighbours nsniinii. Localised spin models are gen¬ 
erally used for insulators. They can be used for metals 
only in the strong coupling limit, when interactions place 
the metal at the verge of a Mott transition. 

Note that magnetism at finite temperatures and many 
experimental features require some three- dimensionality. 
Nevertheless most of these low-energy, localised or multi¬ 
orbital, models are two-dimensional. 


III. ELECTRONIC CORRELATIONS AND (tt, 0) 
COLUMNAR ANTIFERROMAGNETISM 

The nature of magnetism is linked to the strength of 
correlations in the non-magnetic state. Three main pic¬ 
tures have been proposed to explain the (tt, 0) ordering 
of iron pnictides: an instability of the Fermi surface due 
to nesting in a weakly correlated state, antiferromag¬ 
netic exchange between localised electrons, and a double¬ 
exchange like mechanism in a Hund metal state with or¬ 
bital differentiation. In the following we review these 
approaches comparing with experiment and ab-initio cal¬ 
culations. We end the section with a brief discussion. 


A. Weak coupling description of magnetism 

The tendency towards a spin density wave instability 
(SDW) with momentum Q = (tt, 0) linked to the Fermi 
surface shape was noticed EH EH before this order was 
found experimentally |32j I132j . Ab-initio calculations 
predict strong nesting between electron and hole pock¬ 
ets connected by Q, as illustrated in Fig. [^b). Nesting 
between idealised circular pockets leads to a divergence 
in the RPA spin susceptibility at Q. Realistic pockets are 
elliptical and nesting is not perfect, nevertheless a peak 
in the spin response is obtained [271 1133j . 

The antiferromagnetic tendencies emerging from the 
Fermi surface were early related to an s^ supercon¬ 
ducting instability [27]. FRG calculations showed that 
the strong antiferromagnetic correlation also drives su¬ 
perconducting pairing, Fermi surface distortions and or¬ 
bital current order On the basis of RG calculations, 

it was argued that magnetic and pairing instabilities are 
determined by the same interband pair hopping 
When electron and hole pockets are nearly identical the 
SDW instability occurs at a higher T. If the differences 
in shape or size between the Fermi pockets become more 
prominent, the Cooper instability comes first [TTlIT^. 

The suppression of magnetism with doping observed 
experimentally is expected as electron or hole doping 
changes the size of the Fermi pockets and destroys nest¬ 
ing |134j . This is consistent with the non-monotonic de¬ 
pendence on doping of the minimal interaction for the 


onset of magnetism Uonset obtained in Hartree-Fock cal¬ 
culations for a five-orbital model PIITTFt] in which the 
smallest Uonset appears at electronic fillings close to that 
of undoped compounds n ^ 6 , see Fig. |^d). However 
FLEX calculations question the role of nesting as the 
later is suppressed by self-energy corrections m- More¬ 
over Uonset > 1 eV seems too large to claim a nesting 
induced SDW [9^ llQ7j . FeTe, with a similar Fermi sur¬ 
face and believed to be more correlated, shows a different 
magnetic ordering not explicable with nesting arguments, 
see Sec. m 

Changes in the size of the electron and hole pockets 
upon doping modify the momentum at which the RPA 
spin susceptibility peaks and could produce a change in 
the ordering moment. The tendency towards incom¬ 
mensurate ordering is reduced by Umklapp processes 
which enhance the magnetic instability with momen¬ 
tum Q, commensurate with the reciprocal lattice m- 
Nevertheless, signatures of incommensurate spin fluc¬ 
tuations appear in the non-magnetic hole-doped com¬ 
pound [HKl IHH] KFe 2 As 2 and long-range and short- 
range static incommensurate antiferromagnetic order has 
been observed in electron-doped BaFe 2 As 2 [1371 [I38]. 
While the short-range order was initially interpreted as 
an evidence of nesting-induced antiferromagnetism, more 
recent results suggest that it is due to a cluster spin glass 
phase [nmno]. 

In the reconstructed bandstructure of the magnetic 
state small Fermi pockets along the (0, tt) direction have 
been observed in quantum oscillation and ARPES ex¬ 
periments |141I[143] . These pockets are obtained within 
low-energy effective models |144j and in Hartree-Eock cal¬ 
culations of five-orbital models |1Q41110611145] only if the 
interactions are not too large. In this case the topol¬ 
ogy of the bandstructure produces nodes in the SDW 
gap, small Dirac pockets along the T — X direction and 
metallicity m- 

In favour of weak coupling approaches, strong signa¬ 
tures of Hubbard bands have not been observed in X-ray 
experiments P I146IIT48] . ARPES and Quantum oscilla¬ 
tions experiments demonstrate the existence of quasipar¬ 
ticles [28ll3Ql 1112111491 ll5Qj . The shape of the bands is 
reasonably well described by ab-initio methods [HIST]. 
Nevertheless the band mass is considerably enhanced 
m*/m ^ 2 — 3 in undoped pnictides [29] [30]. The low 
ordered magnetic moment, between 0.09jub in NaEeAs 
and 1 /iB in SrEe 2 As 2 [KlIUHl. and the metallicity of the 
magnetic state have been also used in favour of an SDW 
description of magnetism [1041 llQ 6 j . The observed linear 
increase of the spin susceptibility with temperature [3] 
has been explained with Eermi liquid arguments, as a 
consequence of antiferromagnetic interactions [M] or 
due to quasiparticle excitations from a sharp peak in the 
renormalized spectral function located closely below the 
Eermi level [151 j . The spectral features in the magnetic 
state of five-orbital models calculated within Hartree- 
Eock compare well with optical conductivity, photoemis¬ 
sion and neutron experiments only if interactions are 
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weak na uni]. Moderate interactions are also neces¬ 
sary to reproduce the sign of the resistivity anisotropy 
observed experimentally in calculations with band re¬ 
construction as the only source of anisotropy ma , see 
Sec.IVl 

On the other hand, the band reconstruction in NaFeAs 
observed by ARPES involves bands well below the Fermi 
level [86j I152j . Inelastic neutron scattering (INS) and 
X-ray experiments have detected spin waves up to ener¬ 
gies above 200 meV lansMnsj, larger than predicted 
by weak coupling approaches |156j and estimated to¬ 
tal moments larger than those expected in the weak 
coupling model 0 0 UMl EHl USa EH] which do not 
seem to change much with increasing temperatures up to 
T ^ 2T/V, well above the Neel temperature. This mag¬ 
netic moment is not only large in systems with magnetic 
order, but also in compounds which do not order mag¬ 
netically as LiFeAs or Bao. 6 Ko. 4 Fe 2 As 2 |148j . In a Fermi 
surface instability picture only the bands close to the 
Fermi level are expected to be affected by antiferromag¬ 
netism, the spin wave spectral weight is observed only at 
low energies becoming overdamped by electron-hole pair 
excitations at high frequencies |1Q511159j and, except for 
fluctuations effects at temperatures above but close to 
the transition temperature, the magnetic moments are 
formed at Tat. 


B. Antiferromagnetism from localised spins 

An opposite point of view to explain the magnetic 
properties is based on localised spins which interact to 
first and second nearest neighbours with antiferromag¬ 
netic exchanges Ji and J 2 [33l |34l [37[ 11261 I127j . see 
Eq. §. A (7r,0) ground state requires J 2 > Ji/2 as 
antiferromagnetic ordering with momentum Q = (tt, tt) 
is lowest in energy if J 2 < Ji/2. For J 2 > Ti/2 the clas¬ 
sical ground state is selected by J 2 . It consists of two 
interpenetrating square sublattices with Neel order, be¬ 
ing the energy independent on the angle between the two 
sublattices. Columnar order is selected by quantum or 
thermal fluctuations through an order by disorder mech¬ 
anism |M1 [T2H1[T^ . Interactions to first nearest neigh¬ 
bours are frustrated. 

The exchange constants Ji and J 2 derived at large in¬ 
teractions from multi-orbital models for iron supercon¬ 
ductors depend on Hand’s coupling and on the atomic 
orbital configuration laasa, see Sec. |V| The large value 
of J 2 required to fulfill the condition for columnar order 
J 2 > Ji/‘^ is possible due to the role of the As atom in 
the hopping between second neighbours Fe atoms [331 [371. 
The exchange to first neighbours Ji includes both direct 
and indirect (via As) contributions. 

The presence of magnetic moments above and sig¬ 
natures of well defined spin waves throughout the Bril- 
louin zone are commonly used to justify the local moment 
picture [7]. The spin waves measured by neutron scat¬ 
tering can be fitted to a Ji — J 2 Heisenberg Hamiltonian 


but fitting parameters are controversial as (i) they can 
be significantly sensitive to the inclusion of long-range 
couplings |16Qj (ii) the values reported imply Jia » Jib 
with Jib slightly ferromagnetic, while in the tetragonal 
state Jia = Jib- Here Jia and Jib are respectively the 
exchange constants along x and y directions [sa USD. 
The difference between Jia and Jib was justified assum¬ 
ing the emergence of orbital ordering Uyz — Uzxi namely 
different filling of zx and yz orbitals, in the nematic state 
above T/v, and a strong impact of this orbital ordering 
on the exchange constants |161l I162j . Different Jia and 
Jib have been obtained in ab-initio calculations at low 
temperatures [HU [85l [128], see below. A spin localised 
model is generally consistent with orbital fillings not far 
from integer. This is at odds with substantial orbital 
ordering. A different proposal to explain the differ¬ 
ence between Jia and Jib involves a non-negligible value 
of the biquadratic coupling K between nearest neigh¬ 
bour spins. K induces an Ising degree of freedom and 
it is behind the selection of (tt, 0) or (0, tt) order from 
the classical manifold as well as behind the nematic and 
structural transitions above T/v- In a columnar ground 
state the inclusion of K gives Jia = Ji + 2KS^ and 
Jib = Ji — 2KS‘^. For well localised spins, K is gener¬ 
ated by fluctuations [12611127] and it is small. Therefore 
the need of a large value of K to fit the experiments is 
believed to imply deviations from the local moment pic¬ 
ture |13Qj . 

One of the main drawbacks of the localised spin model 
is the size of the ordered magnetic moment (see Table 
in Ref. 0). The crystal field arrangement of iron pnic- 
tides predicts that the 6 electrons may be accommodated 
in five d-orbitals with atomic magnetic moment |5'| = 2 
(high Hand’s coupling) or |5'| = 1 (small Hand’s cou¬ 
pling). These atomic moments are too large compared 
to the magnetic moment measured by neutron diffrac¬ 
tion experiments: y ^ lys in SrFe 2 As 2 and smaller in 
other compounds [5j I155j . 

The early proposal that the frustation of first near¬ 
est neighbour interactions could strongly reduce the pre¬ 
dicted magnetic moment has not been supported by ex¬ 
plicit calculations [T 2 HI UMl UMl MM- The extent of 
frustration depends crucially on the ratios of exchange 
constants and their anisotropies and it is reduced if 
Jia 7^ Jib- The suppression of the magnetic moment 
is small for 1 < |5'| < 3 and it decreases with increas¬ 
ing J 2 |164j . Unrealistically, low atomic spin and | J 2 /T 1 I 
values very close to 1/2 would be required to explain the 
strong reduction of the ordered magnetic moment found 
in experiment. 

The magnetic moments at high temperatures esti¬ 
mated by integrating the spectral weight in neutron scat¬ 
tering experiments is much smaller than expected for lo¬ 
calised electrons prras]. The moments detected by core 
level photoemission spectroscopy (PES) |158j and x-ray 
emission spectroscopy (XES) |157j are larger than those 
measured by INS. However the different results of these 
techniques suggest fast fluctuations of the atomic spin. 
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beyond the localised model |148j . Consistently, the lo¬ 
calised moment approach can neither explain the tem¬ 
perature dependence of the spin excitations |165j nor the 
longitudinal spin excitations detected |166j in BaFe 2 As 2 : 
The spin waves in a localised moment model should be 
purely transverse spin excitations with moments fluctu¬ 
ating perpendicular to the staggered magnetisation. 


C. Insight from ab-initio calculations 

The (tt, 0) order has been widely studied with DFT 
methods. Early calculations resulted in different order¬ 
ings in the ground state: proximity to weak ferromag¬ 
netism [^I167j . (7r,7r) antiferromagnetism |167[I168] and 
(7r,0) stripe antiferromagnetism [32[ [33]. It was then 
clarified that (tt, 0 ) is the magnetic ground state of iron 
pnictides and that pseudopotential approaches have to be 
used with special care, only after being checked against 
all-electron methods [82], because they can fail to iden¬ 
tify the correct ground state due to the small differences 
in energy between the different magnetic orders found in 
these materials. Some calculations were interpreted in 
favour of a weak-coupling origin of magnetism |169j , but 
others arrived to opposite conclusions jSS] [86] . 

The magnetic moment obtained [82] in the (tt, 0) state 
is ^2/i^, much larger than the experimental one 0.35- 
IflB 0 . Such an overestimate of the magnetic moment is 
rare in DFT. A proposal to solve this discrepancy involves 
a high moment ground state with well defined spin den¬ 
sity waves but with magnetic twins and antiphase bound¬ 
aries dynamic on the time scale of the experiment [nni. 

The calculations evidence strong magnetoelastic cou¬ 
pling: (i) The moment is very sensitive to the As-z po¬ 
sition. Small displacements of the As change the Fe- 
moment from 2.0 to 0.5 /ib [331(8^ (8411128] : and (ii) 
the atomic positions, generally well described by GGA 
calculations, show important disagreements with exper¬ 
iments [sa ESI mi]- The discrepancies in the crystal 
structure are smaller in non-magnetic materials while 
spin polarised calculations are in agreement with the ex¬ 
perimental structure in the non-magnetic state [82] [85] . 
Magnetic order has to be included to reproduce the ex¬ 
perimental phonon spectrum [85] 1172] . 

The sensitivity to the approximation used and to tiny 
details of the crystal structure was initially ascribed to 
the itinerancy of the magnetic ground state which was ar¬ 
gued not to be describable in terms of local magnetic mo¬ 
ments [82] . The discrepancies between non-magnetic cal¬ 
culations and the experimental structure were later inter¬ 
preted as an evidence that in experiment there are mag¬ 
netic moments at high temperatures, above Tat, which 
cannot be well captured by DFT methods [85] . 

The value of the exchange constants has been calcu¬ 
lated by DFT using different methods: (i) comparing the 
energies of the magnetic states and mapping them to a 
Heisenberg model, (ii) using the direct spin-flip method, 
appropriate for high-temperatures [85] , and (iii) with lin¬ 


ear response perturbation theory, suitable for low tem¬ 
peratures [85] 112811129] . The J 2 obtained from the mag¬ 
netic state energy is large, as required to stabilize the 
stripe ordering in the localised picture. However, the 
validity of the method was questioned by a perturbative 
calculation which found that different magnetic states are 
stabilised in different atomic orbital configurations [37] . 
see Sec. El 

Direct spin-flip and linear response perturbation the¬ 
ory provide the dependence of the exchange constants on 
the distance. Interactions are short-range in local mod¬ 
els and long-range when magnetism is driven by nesting 
of the Fermi surface. In both approaches, the exchange 
constant along the diagonal direction was found to be an¬ 
tiferromagnetic and very short ranged [85] . Along x and 
y axis, they decay more slowly [85], as 1/R^. Surpris¬ 
ingly, Ji 5 is very different in both methods: it is antifer¬ 
romagnetic and close to Jia at high temperatures (direct 
spin-flip method) and ferromagnetic at low temperatures 
(linear response) [85l 1129] . The later result is compat¬ 
ible with neutron experiments. However, different Jia,b 
at low and high temperatures is at odds with a Heisen¬ 
berg description [85]. The different magnetic moment at 
which DFT calculations converge for different magnetic 
patterns also points against a Heisenberg description of 
magnetism [ 8 ^ . 

LSDA calculations allowed variations on the angle 0 
formed by the magnetic moments on the two Fe inter¬ 
penetrating antiferromagnetically ordered sublattices se¬ 
lected by the classical solution of the Ji — J 2 Heisenberg 
model They found a strong dependence of the en¬ 

ergy on 0 and attributed it to deviations from a Heisen¬ 
berg description which could be accounted for by a large 
biquadratic coupling term K, which, as discussed above, 
can be behind the different values of Jia and J 15 . A 
recent study on several families of Fe-pnictides [174] has 
confirmed that a biquadratic term K is required to fit the 
ab-initio energies with a spin model. K is of the same or¬ 
der of magnitude as the exchange interaction, it depends 
on the Fe-Fe and Fe-As distances and it is substantially 
reduced by hole-doping HU. 

Based on an analysis of the LDA single electron ener¬ 
gies obtained for BaFe 2 As 2 and FeTe with different mag¬ 
netic orderings, Johannes and Mazin argued against the 
binary choice between superexchange and Fermi surface 
nesting in favour of a third mechanism that is neither 
fully localised nor fully itinerant [ 86 ]. They concluded 
that magnetic moments at Fe sites appear due to Hand’s 
rule coupling but that superexchange is not operative. 
The interactions between the Fe moments are consider¬ 
ably long-range. However the role played by Fermi sur¬ 
face nesting in setting the magnetic state is small as for 
both materials the energy gain in the magnetic state is 
largely due to the band reconstruction far from the Fermi 
level [ 86 ] . 
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FIG. 3: (Colour online)(a) Phase diagram of correlations, as measured by the quasiparticle weight as a function of the 
intraorbital interaction U and Hund’s coupling Jh calculated with the U(l) Slave Spin technique for undoped superconductors 
with n = 6 electrons for a five-orbital model [89] proposed for LaFeAsO, adapted from m- Three regions can be differentiated 
according to the value of the quasiparticle weight: (i) A Mott insulating state at large interactions, (ii) a moderately correlated 
metal for small U or Jh with > 0.5, and (iii) a strongly correlated metal with small quasiparticle weight and strong 
orbital differentiation (the Hund-metal), see also Figs. a) and (b). Mott insulators with \S\ = 1 (low Hund) and \S\ = 2 
are separated by a dot-dashed line at large U and small JhIU. (b) Orbital filling versus U calculated with DMFT using full 
rotational invariant spin interaction for the same tight-binding model in (a), Jh/U — 0.25 U < 3 eV and Jh — 3 eV for 
U > 3 eV, adapted from ESI. ( c) Spin fluctuations as a function of U for JhIU — 0.22 for a five-orbital model for FeSe. 
A strong enhancement in the spin fluctuations is observed at the crossover at which the system enters into the Hund metal, 
adapted from [116] . (d) Dependence on doping of the interaction which separates the weakly correlated state from the spin 
freezing state at finite temperature, obtained in DMFT [llOlllll] calculated with the tight-binding model and Jh used in (b). 
A similar doping dependence is expected for Un^nd lower temperatures. With hole-doping Unund approaches the interaction 
Uc at which the Mott transition happens at half-filling. 


D. The role of Hund’s coupling and the orbital 
degree of freedom 

As discussed above, neither weak nor strong coupling 
approaches are able to explain all the experiments. The 
renormalisation of the mass, which in a Fermi liquid de¬ 
scription is the inverse of the quasiparticle weight, is ^ 3, 
emphasizing the importance of correlations. Aside from 
a few exceptions [37[ |75l 117511176j weak and strong cou¬ 
pling approaches ignore the physics emerging from the 
multi-orbital character of iron superconductors. The sim¬ 
plest generalisation of these models allows the coexis¬ 
tence of itinerant weakly correlated orbitals and strongly 
localised ones. However, in the last years it is becoming 
clear that Hund’s coupling plays a key role on the cor¬ 
relations of these materials and other multi-orbital sys¬ 
tems ISa EH ESI Una EHl EHl EZHIHI] • 

1. Correlations induced by Hund’s coupling: Hund metals 

DMFT, Slave Spin and Gutzwiller techniques, able to 
address the renormalisation of the quasiparticle weight, 
have been used to analyze the correlations induced by 
Hund’s coupling in multi-orbital systems [42j |43l I177F 
ll8Qj . Soon after the discovery of superconductivity in 
doped LaFeAsO, LDA+DMFT calculations [39l 1^ flOQ] 
determined this parent compound to be a correlated 
metal far from the Mott transition and not well de¬ 
scribed by either atomic physics or band theory HOD]. 
Hund’s coupling was found to be responsible for the elec¬ 


tronic mass enhancement [39j [42] . At large Jh the cor¬ 
related metal, in the following Hund metal, is built from 
high-spin states with small overlap with single particle 
states liDiHDi rrTDirrsD] . suppressing Z. 

In multi-orbital systems, the Mott transition happens 
not only at half-filling but at all atomic integer fillings. 
The intraorbital interaction Uc for the Mott transition 
depends on the filling. Hund’s coupling modifies Uc [IS 
117^118011182] : (i) At half-filling correlations increase and 
Uc decreases with Jh- (ii) For one electron or one hole 
occupancies Jh reduces correlations, promotes metallic 
behaviour and Uc increases, (iii) Otherwise, Uc is non- 
monotonous with Jh and induces bad metallicity [43l 

diH]. 

The non-monotonic case is relevant to undoped pnic- 
tides with 6 electrons in 5 orbitals. Fig. [^a). With in¬ 
creasing [/, at intermediate and large Jh the system does 
not evolve directly from a weakly correlated metal to a 
Mott insulator |43l IllQl 11161 11181 11771 1179] . Instead, 
there is a crossover between the weakly correlated re¬ 
gion and a strongly correlated Hund metal with small 
Z, see Fig. [^a) an d before the system turns into a 
Mott insulator |43l 111611118] . The interaction Un^nd 
which the Hund metallic behaviour appears depends on 
the number of electrons and on Jh and it can be smaller 
than the bandwidth W [431 11161 1118] . In the Hund 
metal spectral weight is transferred into broad Hubbard 
bands [Tnni[T78l[T83] . 

The reduction of the quasiparticle weight Z originates 
in the suppression of atomic configurations which reduce 
the atomic moment, especially those with double orbital 









































occupancy [43]. It is concomitant with the enhancement 
of spin correlations |40l HU |43l 11161 IlSOj as illustrated in 
Fig-i c), almost saturated in the Hund metal. Hopping 
processes which do not reduce the magnetic moment are 
allowed and the system is metallic 113 nTg. In Fermi 
liquid theory, the quasiparticle weight is the overlap be¬ 
tween the elementary excitations of the interacting and 
the non-interacting systems. The small value of Z in 
the Hund metal is due to the small overlap between the 
locally spin polarised atomic states and single-particle 
states, but it does not necessarily imply strong localisa¬ 
tion [43] of the atomic charge. 

While the spin correlation functions in the Hund metal 
paramagnetic state correspond to a large local moment 
on short time scales, if the electron mobility is high 
enough, this local moment fluctuates very fast so that on 
longer time scales the time-averaged magnetic moment is 
considerably reduced m- With increasing temperature 
the screening time increases and the screened magnetic 
moment at intermediate time scales becomes closer to its 
value at short times |185j . A crossover towards a highly 
incoherent metallic state with local frozen moments takes 
place m- A characteristic temperature, the coherence 
temperature T*, can be defined. Below T* the magnetic 
moments are screened and Fermi liquid behaviour is ob¬ 
served. Above T* the spins freeze |11Q[ 1177] I186j . As a 
consequence, the spin susceptibility evolves from an en¬ 
hanced Pauli susceptibility (typical of Fermi liquids) at 
low temperature to a Curie-Weiss behaviour characteris¬ 
tic of localised spins [42| . A Curie behaviour is expected 
for large screening times |185j . In the frozen moment 
state the self-energy shows a non-Fermi liquid power-law 
dependence [T771 ITMl IW] on frequency and the 
resistivity strongly increases [42|. The coherence tem¬ 
perature T* becomes smaller with increasing Jh |186j . 
Screening large local moments is difficult and T* can be 
very low. 

The value of and of the coherence tempera¬ 

ture T* decrease and the strength of correlations in¬ 
creases as the half-filled Mott insulator is approached uni 
IllOi 11191 11771 I186j . In pnictides this leads to the 
enhancement of the effective mass and the suppres¬ 
sion of the Drude peak Hi nm with hole-doping, see 
Fig. [^d) and Fig. Qc). This has led some authors 
to describe iron superconductors as doped Mott insu¬ 
lators @01 inni nm nzzi, making connection with the 
physics of cuprates. 

Based on DMFT calculations, Ishida and Liebsch 
placed LaOFeAs at the weakly correlated side of this 
crossover mnum]. Slave Spin results situate BaFe 2 As 2 
around the crossover or just above, into the corre¬ 
lated phase |40|. Other LDA+DMFT calculations 
place most iron superconductors in the Hund metal re¬ 
gion [35|. Chalcogenides are believed to be more corre¬ 
lated than pnictides, especially FeTe and alkaline doped 
Ka,Fe 2 -ySe 2 . The later compound shows insulating be¬ 
haviour when vacancies are present [21]. LDA-fDMFT 
describe FeSe as a Hund correlated metal |111[ I188j . 


Based on Gutzwiller calculations it has been argued than 
FeSe and FeTe |116j are Hund metals while using Slave 
Spin it has been proposed that alkaline doped systems 
are close the border of the insulating state m\- 

Experimentally, in most undoped 1111 and 122 FeAs 
compounds the mass enhancement estimated from quan¬ 
tum oscillations and ARPES [2^[3QHTT^I149j is m*/m ^ 
2 — 3. This value is consistent with the quasiparticle 
weight Z at the crossover between the two metallic re¬ 
gions or in the Hund metal one ignn], see Eig. ia). 
Larger mass enhancement were reported for 111 com¬ 
pounds up to m*/m ^ 6 in NaEeAs from ARPES ex¬ 
periments |152j and m*/m ^ 5 in LiEeAs from quantum 
oscillations [nni. EeP phosphides are slightly less corre¬ 
lated [28] with m*/m ^ 1.4 — 2 while for EeSe and EeTe 
the enhancement is between 3.5 — 7 ma. The multi¬ 
orbital character of iron superconductors make optical 
experiments difficult to interpret nevertheless 

the importance of correlations can be also inferred from 
these measurements |192j . 

The characteristic time scale for the screening of the 
magnetic moment in iron pnictides was estimated |185j 
^ 10 fs. This implies that the magnetic moment de¬ 
tected should depend on the experimental technique, be¬ 
ing larger in the case of fast techniques as X-ray, and 
smaller for techniques which average the magnetic mo¬ 
ment in longer time scales such as neutron diffraction, 
/iSR or NMR |148j . In INS experiments the time scale 
is determined by the maximum energy up to which the 
spectral weight is integrated, but it is nevertheless larger 
than in the X-ray techniques PES and XAS. In agree¬ 
ment with this prediction, experimentally the values of 
the magnetic moments estimated by PES measurements 
are much larger than those obtained by neutron scatter¬ 
ing [5l 11481 [T58] . Note, that the comparison between PES 
and neutron measurements is not obvious, as PES deter¬ 
mines the total uncorrelated spin while present inelastic 
neutron scattering measures correlated spin excitations, 
and thus will underestimate the size of the effective spin 
when excitations become diffusive [5]. A temperature 
dependent screening time could be behind the increase 
in the magnetic moment with temperature detected in 
Eei.iTe [l93] 

The spin susceptibility at room temperature is en¬ 
hanced with respect to its bare Pauli spin susceptibility 
value. In 1111 and 122 compounds the magnetic sus¬ 
ceptibility increases with temperature up to, at least, 
T 700 K [3]. This temperature is probably below the 
coherence temperature T* of iron arsenides [1851 I186j . 
The increasing magnetic susceptibility could be a conse¬ 
quence of the antiferromagnetic correlations [3 d fm] 
or due to a sharp peak in the renormalized density of 
states with an effect similar to that of a van Hove singu¬ 
larity On the other hand a crossover is observed in 

the spin susceptibility of EeSe, it increases with temper¬ 
ature below 180 K and it decreases for larger tempera¬ 
tures |194j . The spin susceptibility of EeTeo .92 decreases 
with temperature above the antiferromagnetic transition 
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FIG. 4: (Colour online)(a) and (b) Orbital-dependent quasiparticle weight versus intraorbital interaction U calculated with 
U(l) Slave Spin for Jh jU = 0.10 and Jh jU = 0.25 respectively, with n = 6 electrons as in undoped iron superconductors and 
the five-orbital model [89] used in Fig. [^a). In (b) is strongly suppressed at the crossover between the weakly correlated 
state and the Hund metal. Orbital differentiation is enhanced in the later. Adapted from [118] . (c) Quasiparticle weight Z^^ 
as a function of electronic filling n for the model in (a) calculated with Z 2 Slave Spin technique (f/ = 3 eV and JujU — 0.25). 
Zr^ decreases with hole-doping as half-filling is approached CSj. 


Tat - 70 K [T95]. 

With hole-doping the electronic filling of iron supercon¬ 
ductors approaches half-filling (5 electrons in 5 orbitals). 
Therefore, if undoped iron superconductors are in the 
weakly correlated region a crossover to the Hund metal 
state would be induced by hole-doping. Similarly, if un¬ 
doped iron superconductors happened to be close to the 
n = 6 Mott insulator, they would become Hund metals 
with electron doping. 

Doping dependent correlations consistent with expec¬ 
tations of Hand’s physics have been measured. Differ¬ 
ent techniques evidence that correlations increase when 
BaFe 2 As 2 is doped with holes and decrease when it is 
doped with electrons [40]. Specific heat and spin sus¬ 
ceptibility measurements indicate that correlations in 
KFe 2 As 2 , strongly hole-doped with n ^ 5.5, are stronger 
than in BaFe 2 As 2 on spite of the smaller cation size 
which increases the covalency |117[ 1196] . Large mass en¬ 
hancements have been measured in KFe 2 As 2 by quantum 
oscillations [197] and ARPES [198] . The Drude peak 
in KFe 2 As 2 is strongly suppressed [19911200] when com¬ 
pared with that of BaFe 2 As 2 in agreement with theoret¬ 
ical estimates The specific heat enhancement is 

even larger in CsFe 2 As 2 and RbFe 2 As 2 , isova- 

lent to KFe 2 As 2 . Moreover, half-filled Mn-compounds 
BaMn 2 As 2 and LaMnPO with 5 electrons in the 5 
d-orbitals are insulators [2031 1204] and the strongly 
electron-doped compound BaCo 2 As 2 with 7 electrons in 
the 5 d-orbitals shows a small mass enhacement [205] 
m ^ 1.4. 

Magnetic susceptibility and thermal expansion mea¬ 
surements sni evidence a temperature induced 
coherence-incoherence crossover in KFe 2 As 2 . At low tem¬ 
peratures (T < bOK) the spin susceptibility is nearly 
independent on temperature while it displays Curie- 
Weiss behaviour at high temperatures T > 150iF. Its 
resistivity is very strongly dependent on temperature 
with an extremely large residual resistivity ratios RRR> 
1000 [11711206] . The superconducting transition is rela¬ 


tively broad, suggesting that this large RRR cannot be 
only due to reduced disorder but it is a consequence of the 
temperature induced incoherence uni- Very large RRR 
and non-Fermi liquid behaviour have also been mea¬ 
sured in CsFe 2 As 2 and RbFe 2 As 2 |201[ 1202] . CsFe 2 As 2 
shows paramagnetic behaviour: in all the range of tem¬ 
peratures measured the spin susceptibility increases as 
the temperature decreases [201] . The temperature de¬ 
pendence of all these measurements is consistent with 
KFe 2 As 2 , CsFe 2 As 2 and RbFe 2 As 2 being Hund metals. 
We note that an alternative explanation for the non¬ 
monotonic temperature dependence of the spin suscep¬ 
tibility of KFe 2 As 2 relies on a sharp peak in the density 
of states below the Fermi level [Ml- However, to our 
knowledge, the monotonic increase of the spin suscepti¬ 
bility with decreasing temperature has not been discussed 
within this single-particle picture. 


2. Orbital differentiation in iron supereonduetors 

The orbital-dependent hoppings and the crystal field 
splittings render the five Fe-d orbitals inequivalent with 
consequences on the correlations. LDA calculations pre¬ 
dict for 3z^ — the largest occupation in most pnictides, 
n 3 ^ 2_^2 ^ 1.5 — 1.6, while the other orbitals have fillings 
rij ^ 1.0 — 1.25 [39]. Small differences are found between 
materials [2081 1209] . With increasing interactions the 
population of 3z‘^ — calculated with DMFT decreases, 
that of x‘^ — y‘^ increases inni uni nnn unsi and rixy 
approaches half-filling, see Fig. |^b). In the Mott insu¬ 
lating state at large V, except for very small Jh, 
is doubly occupied while the other orbitals are half-filled. 

Hand’s coupling enhances the differences in filling and 
quasiparticle weight among the orbitals and can lead to 
orbital selective Mott transitions (OSMT) [39f[4T] . In 
the state with moderate correlations in Fig. (a) the 
differences among the orbital dependent Z^ are small, 
see Fig. I^a). However, according to DMFT, Gutzwiller 
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and Slave Spin [40l IllQi lllli 11161 11181 I189j . they are 
substantial in the correlated Hund metal, see Fig. [^b). 
The t 2 g orbitals, specially xy, become more correlated 
than the eg orbitals [M14T1 iTTom2i [ml 
I119[I189[ 121011211j . Orbital differentiation increases with 
hole-doping [40l I119j and could be behind the sensitivity 
to structural details maiiis]. 

Among the experimental signatures of orbital differ¬ 
entiation are: (i) ARPES and quantum oscillations ex¬ 
periments have reported important differences in or¬ 
bital correlations in KFe 2 As 2 |197[ 1198] : (ii) The de¬ 
viations from LDA predictions of the Fermi surface of 
LiFeAs, LaFePO and LiFeP observed experimentally 
have been adscribed to the orbital dependent mass en¬ 
hancement and the interaction dependent orbital occu¬ 
pancy |212[|213| : (hi) ARPES experiments in 122 chalco- 
genides have been interpreted in terms of a temperature 
induced OSMT [dTlfTS^. 

The spin excitations measured in neutron experi¬ 
ments [71 1136] seem consistent with the presence of both 
itinerant and localised electrons. Itinerant electrons, 
sensitive to changes in the Fermi surface with doping, 
are assumed to be responsible for the low-energy excita¬ 
tions while localised electrons account for the high-energy 
spectral weight. The integral of the spectral weight mea¬ 
sured with INS, which gives the local ffuctuating moment, 
is however strongly reduced in the strongly hole-doped 
KFe 2 As 2 , with respect to undoped BaFe 2 As 2 |136] . This 
result is not easy to understand within a Hand’s metal 
description as a larger atomic moment is expected when 
half-filling is approached. We are not aware of X-ray 
estimates of the magnetic moment in this compound. 

The orbital differentiation has consequences on the in¬ 
terpretation of experiments. In particular the suppres¬ 
sion of the Drude peak spectral weight in optical con¬ 
ductivity experiments and that of the kinetic energy are 
not equal, even in a Fermi liquid description mg. Note 
that the identification of the Drude peak spectral weight 
is also affected by interband transitions |119[ 1191] . 

The dependence of interactions on frequency |186] . 
non-local correlations beyond single site approxima¬ 
tions mu, and deviations from the nominal filling not 
captured by five-orbital models [1881 1215] have been 
claimed to have an effect on the correlations in iron super¬ 
conductors, but they are beyond the scope of this paper. 


3. The role of Hund’s eoupling and orbital differentiation 
in the magnetie state of iron pnietides 

Hand’s coupling and the orbital degree of freedom have 
a very significant effect on the magnetic state of iron su¬ 
perconductors. Fig. 1^ (a) shows the magnetic phase dia¬ 
gram as a function of the interaction, calculated within k- 
space Hartree-Fock and enforcing (tt, 0) ordering [9^1107] 
(see Sec|V|for a discussion of other magnetic states). 

NM, in grey, is the non-magnetic state. The phase in 
red is a low moment (LM) state, with antiparallel mag¬ 


netic moments, which violates Hund rule. In the blue 
region, all the magnetic moments point in the same di¬ 
rection and Hand’s rule is satisfied. At the right of the 
solid line the system is insulator, metallic at the left. In¬ 
sulating phases with S = 1 and S' = 0, not described 
here, are found at t/ > 4 eV and JhIU < 0.02 [92] . 

The LM state, in red, is characterised by a small mag¬ 
netic moment, see Fig. [^a), which originates in the par¬ 
tial cancellation of the magnetic moments of different or¬ 
bitals [92] . This state arises from the anisotropy of hop¬ 
ping parameters and magnetic exchanges between orbital 
magnetic moments. It emerges as a way to eliminate the 
frustration in the exchange interactions of first nearest 
neighbours [92]. This phase was initially proposed to 
explain the low magnetic moment in LaFeAsO [92]. In 
Hartree-Fock it appears for Jh values smaller than ex¬ 
pected in real materials [92| |98l [TnH] . With increasing 
Jh , its stability is reduced as interatomic and intraatomic 
exchanges compete. A low magnetic moment state vio¬ 
lating Hand’s rule was also obtained in LDA+U calcu¬ 
lations for different values of Hand’s coupling |2I6[ 1217] . 
A Gutzwiller calculation, however, restricted this phase 
to a small region of interaction parameters mg. 

Although the values of the interactions in Fig.|^a) and 
Fig. [^a) are not comparable (the effect of interactions is 
overestimated by Hartree Fock and underestimated by 
Slave Spin), there are important similarities between the 
two phase diagrams. The boundary of the insulating blue 
state in Fig.|^a) shows a non-monotonic dependence on 
interactions similar to that of the Mott insulator in Fig. 
I^a). In the insulating blue phase in Fig. [^a), at the 
right of the solid line, a gap is open at the Fermi level in 
all the orbitals. x‘^ — is completely filled and has zero 
magnetic moment while the other orbitals are half-filled 
and completely spin polarised. These fillings and orbital 
polarisations are also found in the Mott insulator in Fig. 
I^a). Moreover, as discussed in Sec. |v|the competition 
with (7r,7r) ordering found in the insulating state com¬ 
pares very well with expectations based on a mapping to 
a Heisenberg model. Hartree-Fock is a mean field tech¬ 
nique and cannot address the localisation of electrons, 
however these similarities suggest a connection between 
the insulating states Fig. |^a) and Fig. |^a). 

Similarly to the non-magnetic state, the magnetic 
phase diagram does not evolve trivially with increas¬ 
ing U. Changes in the orbital filling are observed, in 
Fig. ^b). 3z‘^ — is partially emptied, xy approaches 
half-filling and the filling x^ — y‘^ increases. This ten¬ 
dency is consistent with the one observed in the non¬ 
magnetic state, see Fig. ib). The novelty in the (tt, 0) 
magnetic state is that the tetragonal symmetry break¬ 
ing induces orbital ordering. The filling of the zx orbital 
increases and differs from that of yz^ which approaches 
half-filling. Half-filled xy and yz show a gap at the Fermi 
level, while partially-filled zx^ x^ — y^ and 3z‘^ — do 
not, see Ref. ung. 

As U is raised, the magnetic moment increases first 
slowly and then faster at see Fig |^c) and 
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FIG. 5: (Colour online)(a) Hartree-Fock (tt, 0) magnetic phase diagram as a function of intraorbital interaction U and Hund’s 
coupling Jh for undoped systems with n = 6 electrons using the hve-orbital model proposed in |31j (slightly different to that 
used in Fig. [^, see text, (b) Orbital Filing in the Hartree-Fock magnetic state as a function of U for Jh/U — 0.25. The 
dot-dashed line marks U^ag^ see Bna. As in the DMFT calculations in Fig. |^b), with interactions 3z^ — is emptied, the 
Filing of increases and xy approaches half-Flling. The novelty in the magnetic state is the different Filing of zx and yz 

consequence of the tetragonal symmetry breaking, yz approaches half-Flling and Uzx increases, (c) Magnetic moment in the 
Hartree-Fock (tt, 0) antiferromagnetic state as a function of U for JujU — 0.25. It shows a strong enhancement at U^ag which 
resembles the enhancement of spin Fuctuations in the non-magnetic state, see Fig. [^c). (d) In black U^ag^ fhe crossover line 
which separates the weakly correlated magnetic state and the one with strong orbital differentiation, and in red Uonsets the 
interaction at which (tt, 0) magnetism sets in, calculated within Hartree-Fock approximation as a function of electronic Filing 
for Juju — 0.25. Uonset is non-monotonous with doping and minimum for n ~ 6 as expected in a nesting scenario. On the 
contrary, U^ag decreases with hole-doping, similarly to U*. All Fgures after HUT]. 


Ref. |lQ7j . This boost coincides with the establishement 
of strong orbital differentiation leading to two distinct 
regions in the magnetic phase diagram HU- At the left 
of the crossover dashed line in Fig. I^a), the mag¬ 

netic moment is small and all the orbitals oehave as itin¬ 
erant, none of them is gapped. At the right of the 

magnetic moment is large and there is strong orbital dif¬ 
ferentiation: zx^ 3z‘^ — and x‘^ —y‘^ are itinerant while 
xy and yz are half-Flled and gapped m\- 

Strong dependence of the magnetic moment on inter¬ 
actions has also been found in Gntzwiller |218j , and Vari¬ 
ational Monte Carlo sm studies of Fve-orbital mod¬ 
els. The enhancement of the magnetic moment at U^^g 
is reminiscent of the spin fluctuations at Both 

crossover interactions and U^^g have similar de¬ 

pendence on [/, Jh and electronic Filing. They decrease 
with increasing Jh and with hole-doping while they in¬ 
crease with electron doping. This behaviour is illustrated 
for Fxed JnjU in Figs. (a) and (d) and|^a) and (d). 
We interpret that in the magnetic state U^^g is pick¬ 
ing up the underlying physics which emerges in the non¬ 
magnetic state at As discussed above, the values 

of and U^^g in both Fgures should not be com¬ 

pared as they have been obtained using different tech¬ 
niques. 

Note in Fig.j^d) the qualitatively different dependence 
on doping of U^^g and Uonsets the minimal interaction re¬ 
quired for the onset of magnetism. While U^^g decreases 
monotonically as the electronic Filing is reduced, Uonset 
is non-monotonic with doping and shows a minimum 
around n ^ 6. The doping dependence of Uonset and the 
connection between U^^g and suggests a nesting¬ 

like origin of magnetism, similar to the one described in 


Sec. Ill A for Uonset < U < U^^g. Interestingly, this 
nesting-like region is not present in the magnetic phase 
diagram calculated with the tight-binding corresponding 
to a squashed FeAs tetrahedra, with a Fermi surface char¬ 
acterised by an enhanced ellipticity of the electron pock¬ 
ets and lack of nesting |1Q711219j . 

The orbital differentiated region in the (tt, 0) Hartree- 
Fock phase diagram for U > U^^g seems connected with 
the orbital differentiated Hund metal in Fig. Hartree- 
Fock cannot address the renormalisation of the quasipar¬ 
ticle weight, nevertheless the observed orbital differen¬ 
tiation suggests a description of magnetism in terms of 
itinerant (zx, 3z^ — r^, x‘^ —y‘^) and localised {xy and yz) 
orbitals. This is consistent with 3z^—r^ and x^—y^ being 
the less correlated and xy the most correlated orbitals in 
DMFT and Slave Spin calculations in the non-magnetic 
state, see Sec. |IIID 2| In the non-magnetic state, zx and 
yz are degenerate and show intermediate correlations. 
The results in the magnetic state suggests that the sym¬ 
metry breaking and difference in Filing between these 
two orbitals make yz more correlated and zx more itin¬ 
erant [Ml- To our knowledge the orbital differentiation 
between zx and yz in the magnetic state has not been 
addressed so far with techniques beyond Hartree-Fock. 


The nature of the magnetism in the orbitally differenti¬ 
ated region differs from the weak and strong coupling ap¬ 
proaches discussed before. A double-exchange-like model 
was introduced |lQ7j to explain the stability of the (tt, 0) 
ordering. Localised xy and yz favour an antiferromag¬ 
netic state with (tt, 0) or (0, tt) momentum depending on 
parameters, see [Ml. The stability of (tt, 0) magnetism 
cannot be explained only considering the localised or¬ 
bitals, see [Ml for a discussion. Itinerant zx, x^ — y^ 
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and 3z^ — r^, with 4 electrons, promote a ferromagnetic 
state to gain kinetic energy, specially in the y direction 
with larger hopping integrals [3TJ [731 EH • Therefore 
(tt, 0) ordering arises as a compromise between the anti¬ 
ferromagnetic tendencies of the localised orbitals and the 
ferromagnetic tendencies, larger in the y direction, of the 
itinerant ones HOT] zx^ — y^, and 3z^ — 


A crossover is also observed in the spin wave spectrum 
of the Hartree-Fock magnetic state calculated including 
all RPA bubbles and ladder diagrams: the spin excita¬ 
tions evolve from broad low-energy modes at weak in¬ 
teractions to sharply dispersing spin waves prevailing to 
higher energies at interaction strengths at which yz and 
xy are half-filled |159j . These spin waves, present up to 
high-energies, are reminiscent of those observed in some 
neutron experiments. 


As discussed in Sec. Ill C , ab-initio calculations were 
interpreted in terms of magnetic moments at the Fe 
sites which appear due to Hand’s rule coupling [86], are 
present in the non-magnetic state [85], and show long 
range interactions on spite of the small role played by 
Fermi surface nesting in setting the magnetic state [86] . 
This proposal is compatible with the system being in 
the Hund metal state. Ab-initio calculations are mean- 
field techniques similar to Hartree-Fock. A possible in¬ 
terpretation of the large magnetic moments ^ 2 ob¬ 
tained in ab-initio calculations and its sensitivity to the 
As position could be that iron pnictides are in the orbital 
differentiation region and the line is crossed with 

small changes in the As position. Monte Carlo calcula¬ 
tions with ab-initio predicted interaction values, in fact, 
place different materials at different sides of the U^^g 
crossover uni. Undoped and electron-doped 122 mate¬ 
rials with larger resistivity in the ferromagnetic direction 
are not expected to be very deep in the orbital differ¬ 
entiation region as orbital-itinerancy along this direction 
would favour an opposite resistivity anisotropy 1 145] . see 
Sec. [El 

Following the discussion above, in the Hund metal the 
local moment fluctuates very fast and the time-average 
moment is reduced The screened moment is 

the one which can be magnetically ordered at low tem¬ 
peratures. Temporal spin fluctuations not captured by 
mean-field methods like Hartree-Fock and DFT would 
explain the reduced size of the magnetic ordered moment. 
This is consistent with X-ray measurements which show 
magnetic moments equal to 2.1//^ in SrFe 2 As 2 and 1.3ys 
in CeFeAsO, respectively larger than their ordered mo¬ 
ments IfiB and 0.8/rB, and than those estimated by INS 
in other pnictides [5l 1148] 1158] . 

The spectroscopic properties of the Hartree-Fock mag¬ 
netic state, such as optical conductivity or Raman scat¬ 
tering |1Q4[ 1145] 1219] , are mostly sensitive to the size of 
the magnetic gap which is directly related to the mag¬ 
netic moment. Therefore the spectrum in the magnetic 
orbital differentiation region could be significantly differ¬ 
ent to Hartree-Fock predictions if temporal fluctuations 
are included in the calculations and such modifications 


could alter the conclusions regarding the region of pa¬ 
rameters compatible with experimental results. 

In the double-exchange model the value of Tjv, its in¬ 
teraction and doping dependences are unknown, n = 6, 
the filling of undoped pnictides at which T/v is maxi¬ 
mum in 122-compounds, does not play any special role 
in the Hand’s metal description. It is not clear whether 
T/v should follow the same dependence on doping as the 
strength of correlations as with hole-doping zx is emp¬ 
tied losing its itinerant character. Moreover, in single 
band models, Tn depends non-monotonically on interac¬ 
tions. TTV can be also influenced by the competition with 
other magnetic orderings, in particular with (i) (7r,7r), 
observed in the chromium-doped compound and ener¬ 
getically favourable with hole-doping, or with (ii) double 
stripe or staggered-dimer orderings |37l 110811220] , found 
in models for iron pnictides and chalcogenides, see Sec.[V| 
Studying Tat as a function of interactions and doping, 
with techniques beyond mean-field, would be very useful. 
This analysis should, nevertheless, use three-dimensional 
models. 


E. Discussion 

Previous discussion leaves us without arguments to 
approach (7r,0) magnetism and correlations in iron ar¬ 
senides by considering only local moments: The materi¬ 
als are metallic, techniques which sample the system in 
different time scales report different magnitudes of the 
magnetic moment, and the ordered moments are small 
and cannot be justified on the basis of frustration. More¬ 
over, longitudinal spin excitations have been detected, 
and the parameters used to fit the spin waves involve 
strong anisotropies, inconsistent with the almost com¬ 
mensurate fillings expected in a localised approach, or a 
large biquadratic term typical of systems which deviate 
from localisation, as also found in ab-initio calculations. 

On the other hand, there are clear signatures of strong 
correlations and of the role played by Hand’s coupling in 
iron pnictides. The electronic bands are strongly renor¬ 
malized. The mass enhancements reported are compara¬ 
ble to those found in the Hund metal state in Fig. |4](b) 
and increase with hole doping. Orbital differentiation, 
i.e. differences in the correlation strength of different 
orbitals, is observed and becomes more pronounced in 
hole-doped systems. There is evidence of magnetic mo¬ 
ments far above the Neel temperature. Their value is 
larger than the one measured by neutron diffraction in 
the ordered state. These moments are present not only 
in systems which order magnetically at low temperatures 
but also in systems which do not order. The presence of 
magnetic moments at high temperatures is in agreement 
with the need of including magnetism to reproduce the 
crystal structure and the phonon spectrum. 

Deep in the Hund metal region a Curie-like tempera¬ 
ture spin susceptibility originating from a small coherent 
temperature T* and larger conductivity in the ferromag- 
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netic direction are expected. Such dependence is not ob¬ 
served in the undoped systems but it has been found in 
hole-doped CsFe 2 As 2 in the whole range of temperature 
measured and in KFe 2 As 2 above 150 K. The anisotropy 
of the resistivity agrees with the experimental one only 
for moderate interactions, see below. 

Therefore, our conclusion is that in undoped FeAs sys¬ 
tems there is evidence of the formation of atomic mo¬ 
ments due to Hund coupling, but the systems are not 
deep in the Hund metal state, i.e. undoped arsenides 
should be placed around the crossover between 

the Hund metal and the weak coupling regions. Hole¬ 
doping moves them deeper into the Hund metal state. 
This conclusion is consistent with theoretical estimates 
which place undoped iron pnictides around the crossover 
U^und non-magnetic state [4Q l [TTOl [TT^ [T86] 

and around in the magnetic state |114j . It is also 

consistent with the strong sensitivity of the magnetic mo¬ 
ment to the As position found in ab-initio calculations. 
Note that while this crossover is relatively sharp in the 
slave spin calculations plotted in Fig. and computed 
including only density-density interactions IMIIIHIIIIH], 
DMFT calculations which treat the rotationally invari¬ 
ant Hand’s interaction find a smoother crossover, see 
Ref. pro] . 

Experimental evidence and theoretical calculations 
show that chalcogenides are more correlated. In partic¬ 
ular, the spin susceptibility of FeSe and FeTe decreases 
with increasing temperature, respectively above 180 K 
and 70 K. Mass enhancements are also larger and an 
anomalous temperature dependence of the magnetic mo¬ 
ment was found in FeTei.i. Moreover the magnetic order 
in FeTe differs from that predicted by Fermi surface ar¬ 
guments. In our view, 11 compounds can be classified as 
Hund metals. The high-moment found in alkaline-doped 
chalcogenides is expected only deep in the Hund metal 
region or very close to the Mott insulating state. 

Finally we note that the origin of (tt, 0) magnetism is 
different at both sides of U^^g\ nesting-like for smaller 
U and double-exchange for larger U. To our knowledge, 
there are no studies of the interplay between these two 
mechanisms at the boundary U^^g. Whether one of them 
dominates, they cooperate or compete is not clear to us. 


IV. ANISOTROPY AND THE 
SPIN-ORBITAL-LATTICE ENTANGLEMENT 

Besides the discussion on the role of correlations in 
the columnar antiferromagnetic ordering, another ongo¬ 
ing debate is the connection of this ordering with the 
nematic state. Q = (tt, 0) columnar ordering is linked 
to the nematic state since it breaks the C 4 symmetry 
of the lattice down to C 2 symmetry. This symmetry is 
also broken by the structural transition and by orbital 
ordering. Consequently, nematicity could have an or¬ 
bital, spin or lattice origin and there is a strong debate 
about its mechanism [ 22111222 ]. The difficulty in decid¬ 


ing between scenarios stems in that they all break the 
same symmetry, which implies that as soon as one sym¬ 
metry is broken in one channel it will be also broken 
in the other two channels. Different experiments seem 
to indicate that nematicity is electronic in origin gsi- 
Ha EH ES ESI EH leaving as candidates the orbital and 
spin degrees of freedom. In the orbital driven nematicity 
scenario [asiizalizi orbital fluctuations provide the pair¬ 
ing mechanism for superconductivity in contrast to the 
spin fluctuation pairing mechanism for the spin driven 
nematicity iniiia. 

Experimentally, the resistivity is smaller for the 
antiferromagnetic direction than in the ferromagnetic 
direction |44f[48] . Counterintuitively, this resistivity 
anisotropy grows with doping and is largest close to the 
border of the antiferromagnetic phase when the mag¬ 
netism and the structural distortion are weaker. The 
opposite resistivity anisotropy has been found in a hole- 
doped material [63]. Anisotropy has also been measured 
in optical conductivity ESMl EH I223j . elastic shear 
modulus @3 EO], Raman m, scanning tunnelling mi¬ 
croscopy [ 52 H 54 ] . magnetic torque [55] and inelastic neu¬ 
tron scattering experiments [5614581 l65] . Signatures of 
orbital-dependent reconstruction of electronic structure 
in the magnetic and non-magnetic states are found in 
ARPES experiments [66ll68l 122411225j and in X-ray ab- 
sortion spectroscopy [69] . 

In addition, there is no agreement about the rel¬ 
evance of impurity scattering in the x/y resistivity 
anisotropy [48l |53| |58l 17611791 1226] . Upon annealing, 
Ishida et al. m find that the resistivity anisotropy of 
BaEe 2 As 2 nearly vanishes while significant anisotropy re¬ 
mains in Co-doped compounds. Similar results are ob¬ 
tained in Ba(Eei_a,Rua,) 2 As 2 [226] . In accordance with 
these experiments scanning tunneling microscopy [52]- 
[54] report anisotropic spatial extended defects, called 
nematogens, at low temperature. However, Kuo and 
Eisher in Ref. [48] find that the strain induced resis¬ 
tivity anisotropy in BaEe 2 As 2 , Ba(Eei_a,Coa,) 2 As 2 and 
Ba(Eei_a,NU) 2 As 2 is independent of disorder in the 
tetragonal phase and it is appreciable in the magnetic 
state of BaEe 2 As 2 . Eurthermore, the nematic suscepti¬ 
bility shows clean critical exponents m, suggesting ir¬ 
relevance of disorder. Theoretically, nematogens around 
impurities have been found in models of iron pnictides 
driven by orbital [223 or spin correlations [228] . On 
the other hand, in calculations in the magnetic phase 
of multi-orbital models the role of the topology and mor¬ 
phology of the Eermi surface in the resistivity anisotropy 
has been emphasised without and with impuri¬ 

ties [229] . 

Hints about the strong interrelation between the spin 
and orbital degrees of freedom come from first principle 
calculations and mean field approaches in multi-orbital 
Hamiltonians. It has been found that orbital ordering is 
present when the (tt, 0) AE sets in in realistic models for 
iron pnictides [TS] 11131II451123011232] . However, to 
the best of our knowledge, the spontaneous generation of 
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FIG. 6: (Color online) (a) U vs J/U magnetic phase diagram superposed to the magnetisation showing the non-magnetic 
(black) and magnetic phases (LM) and (HM). LM stands for low magnetic moment with violation of Hand’s rule and HM 
stands for parallel magnetic moments. The magnetic moment reaches values as high as 3.4 fiB- The color scale emphasizes 
the region with magnetic moments smaller than 2 {ib- (b) Drude ratio Dx/Dy. The experimental anisotropy corresponds to 
Dx/Dy > 1. Largest values of DxjDy > 1 occur at small values of the magnetisation, (c) Orbital ordering riyz — rizx- The 
orbital ordering anticorrelates with Dx/Dy > 1. The figure is reproduced from [145] . 


orbital ordering has not been found without magnetism 
in these type of models. To study the consequences of 
orbital ordering in the nematic phase a phenomenological 
crystal field term breaking yz/zx degeneracy has to be 
included in the Hamiltonian (see for example [36l [79l 
[23311234]). 

Analysis of the Drude weight in multi-orbital mod¬ 
els has helped clarify the mechanism of the resistivity 
anisotropy. In [145] the ratio of the Drude weights 
Dx/Dy is calculated in the magnetic state assuming an 
isotropic scattering rate. Fig. [^a) compares the mag¬ 
netic moment (a), the Dx/Dy ratio (b) and the orbital 
order Uyz — Uzx (c) in the U vs Jh phase space. LM 
and HM in Fig. [^a) stand for the phases which violate 
and fulfill Hand’s coupling, respectively, see the red and 
blue regions in Fig. The black color corresponds to 
the non-magnetic region. In Fig. [^b) a wide range of 
values (0.2 < Dx/Dy < 1.7) is found depending on the 
value of the interacting parameters. The experimentally 
observed transport anisotropy for undoped and electron- 
doped compounds [444(48] corresponds to Dx/Dy > 1. 
These values coincide with the smallest values of the mag¬ 
netic moment and appear within a great part of the LM 
and in the HM phases restricted to the nesting-like re¬ 
gion and to the orbital differentiated region close to the 
crossover to the nesting-like in Fig. This is believed 
to be the physical region for undoped iron arsenides, see 
Sec. jlllDj The result is in accordance with experiments 
in undoped and electron-doped samples since the resistiv¬ 
ity anisotropy is largest close to the magnetic transition 
with small magnetic moment. 

DjDy > 1 appears when the orbital ordering is small 
{Wyz — '^zx\ < 0.05) and Dx/Dy decreases with increas¬ 
ing orbital ordering in most part of the orbital differen- 
ciated region of the phase diagram. Thus, within this 
approximation, orbital ordering cannot be the origing 
of the transport anisotropy since it seems to favour the 


opposite anisotropy. This result is in agreement with 
other works [23511236] . Deep in the strong orbital dif¬ 
ferentiation region, Dx/Dy < 1 is related to the itin¬ 
erancy of the zx^ x‘^ — and 3z‘^ — orbitals in the 
^-directi on whi ch provides stability to the (tt, 0) AF state 
(see Sec. HID) [107] . 


Small orbital ordering and magnetic moments in the 
HM region with Dx/Dy > 1 are consistent with recon¬ 
structed Fermi surfaces inferred from the Q = (tt, 0) dou¬ 
bling of the antiferromagnetic unit cell. In Fig. the 
reconstructed bands in the antiferromagnetic phase are 
shown at low energy. The magnetic bands are calcu¬ 
lated for U = 1.6 eV and J = 0.254/ corresponding to a 
magnetic moment of ^ 0.9//^ appropriate for the exper¬ 
imental values of 122 compounds. An orbital splitting of 
the zx/yz weight at the F pocket can be appreciated as 
observed in ARPES |66] . The connection of the F pocket 
with the X pockets by the Q = (7r,0) antiferromag¬ 
netic vector gives rise to the metallic nodal spin density 
wave [144] with the formation of two Dirac cones along 
V ^ X. These Dirac pockets govern the anisotropy 
in the yellow region of Fig. [^b). In the orange region 
with the highest anisotropy additional small hole pockets 
appear around the Y point of the Fermi surface (Lifsthiz 
transition). The particular form in which the Lifsthiz 
transition arises depends sensitively on the low energy 
physics of each tight-binding model but it will never¬ 
theless influence the Drude anisotropy. This situation 
can be inferred from Fig. shifting upwards the upper 
folded band with zx component close to Y. The ratio 
value Dx/Dy also depends on the shape of the Fermi 
surface. Therefore in this scenario, so-called ’’Fermi sur¬ 
face anisotropy”, the anisotropy is strongly influenced 
by the topology and morphology of the Fermi surface. 
Other multi-orbital Hubbard models [23711238] and effec¬ 
tive spin-fermion models [23911240] for appropriate values 
of the interactions and ab-initio calculations im, have 
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FIG. 7: (Color online) Energy bands in the (tt, 0) antifer¬ 
romagnetic state close to the Fermi level for U = 1.6 eV 
and Jh = 0.251/. Linewidths and colors reflect the orbital 
content yz=Ted, zx=green, xy=h\ue, 3z^ — r^=orange and 
— y^=black. Note that not all the orbital content is visible 
due to overlapping curves. The optical (axx, cFyy) and Raman 
transitions {xBig^ XB 2 g) shown. The figure is taken from 
Ref. HH]. 


also found that the Fermi surface anisotropy corresponds 
to the experimental anisotropy. 

In Ref. |229j the in-plane resistivity and the Drude 
weight were calculated in the antiferromagnetic state 
adding non-magnetic impurities to a five-orbital model. 
It was concluded that the anisotropy is enhanced by the 
effect of the impurity scattering. The scattering rate 
is strongly influenced by the change of the topology of 
the Fermi suface while the Drude weight is influenced 
by the topology and morphology of the Fermi surface 
in accordance to Ref. |145j . Lifsthiz transitions driven 
by temperature in the reconstructed antiferromagnetic 
Fermi surface have also been claimed to be responsible 
for the anisotropy in the resistivity in the impurity sce¬ 
nario m- In summary, both scenarios, disorder driven 
anisotropy and Fermi surface anisotropy, are probably 
relevant to understand the anisotropy in the magnetic 
state and strongly depend on the topology of the recon¬ 
structed Fermi surface. 

The Drude anisotropy extends to finite frequencies as 
seen in optical conductivity experiments j59f[64l 122 ^ . 
Optical conductivity along the antiferromagnetic (Txx{^) 
and ferromagnetic (Jyy{oo) directions show different inten¬ 
sity and peak frequencies. The Raman spectrum in the 
magnetic state shows signatures and peaks at energies 
similar to those found in optical conductivity [24111242] . 
Calculations in the Hartree-Fock approach for magnetism 
show that the different frequencies at which and (jyy 
peak reflect the magnetic gaps at Y and X electron pock¬ 
ets respectively. On the other hand, both interband tran¬ 
sitions are active for the Big polarisation of Raman ex¬ 
periments (see Fig.[^. These transitions are based in the 
velocity and Raman vertices studied in Ref. m and de¬ 
pend on the symmetry of the orbitals involved in the in¬ 
terband transitions. The dependence on /c-space reflects 
the underlying lattice |219j . This observation can help 
interpret optical and Raman data in the magnetic state. 


In the nematic state the situation is even more complex 
since going beyond mean field like descriptions in five- 
orbital models or ab-initio calculations is very computa¬ 
tionally demanding. Two reviews of the nematic state by 
Y. Gallais and I. Paul and by A. Boehmer, J. Schmalian 
and C. Meingast are presented in this volume [24311244] . 
In this context Landau approaches have been shown to 
be very useful to understand the interplay between the 
structural and magnetic degrees of freedom and to cal¬ 
culate the response functions 

Within the spin-nematic scenario [63l |79l 1233] . the re¬ 
sistivity anisotropy originates in the scattering rate af¬ 
fected by the anisotropic spin fluctuations. In the weak 
coupling limit the spin nematic phase is derived from a 
microscopic Hamiltonian with hole and electron pockets 
without structure in the orbital degree of freedom [233] . 
The results crucially depend on the ellipticity of the 
electron-pockets vanishing for circular electron Fermi sur¬ 
faces [23311246] . However, these Fermi pockets posses a 
non trivial orbital composition as shown in Fig. The 
particular arrangement of the yz and zx orbitals arises 
because under a 90 degree rotation the two orbitals trans¬ 
form as \xz) \yz) and \yz) —\xz) [250] . Con¬ 

sequence of this singular C 4 symmetry is a non triv¬ 
ial topology manifested in the vorticity two of the F 
pocket [HlElOlESI]. 

Recently an effective action has been proposed for the 
magnetic instability derived from a multi-orbital Hamil¬ 
tonian m- The Landau coefficients depend on the or¬ 
bital content, Hubbard and Hund’s coupling. In this ap¬ 
proach it is uncovered that the orbital degree of freedom 
changes the spin-nematic scenario in an essential way due 
to the connection between the F pocket with vorticity 2 
and the yzjzx composition of the XjY electron pockets. 
It is found that in the nematic state the yzjzx orbital 
degeneracy is lifted due to spin fluctuations without the 
need of adding a phenomenological crystal field term as 
in previous studies [31 [ 73 1133]). These results extend 
the mean field findings of orbital ordering in the mag¬ 
netic phase to the nematic state [7311^ irm 114^ [23^ 
1232] . It is also found that the spin fluctuations in the ne¬ 
matic phase present anisotropic momentum dependence 
X(7r,0)(q) ~ X(7r,0)(0) + + ay(^y without invoking el- 

lipticity of the electron pockets. This anisotropy is nec¬ 
essary to understand neutron scattering experiments at 
low energy [65l I252H255] and is consistent with interpre¬ 
tations using ab-initio calculations [25311255] . 

The (^4 symmetry breaking in the (tt, 0) antiferromag¬ 
netic state is also reflected in the anomalous behavior 
of specific phonons [ 6 OI 12411125611266] . Particularly in¬ 
teresting is the case of the Aig As-phonon resonance 
in the Raman spectroscopy upon entering the magnetic 
phase. The Aig As-phonon involves arsenide atoms vi¬ 
brations along the 2 )-axis. In the tetragonal phase the 
Aig As-phonon is active neither in Big nor in B 2 g. Big 
and B 2 g polarisations are sensitive to anisotropic re¬ 
sponses along the Fe-Fe bond and Fe-As bond respec¬ 
tively. A strong signal, which cannot be explained by 
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the orthorombic distortion alone, appears just in Big Ra¬ 
man response when undergoing the magneto- 

structural transition . In Ref. |268j the Raman response 
is evaluated in the non-magnetic and in the magnetic 
states. It is found that the magnetic phase is able to in¬ 
duce a Raman intensity in the Big Raman polarisation 
due to the C 4 symmetry breaking but the intensity is very 
sensitive to the way in which the coupling between the 
electrons and the phonons is introduced. With similar 
symmetry arguments, a finite Te-phonon intensity in the 
B 2 g polarisation was predicted to appear in the double 
stripe magnetic state of FeTe |268j . 


V. MAGNETIC SOFTNESS 

Although the (tt, 0) columnar magnetic ordering is the 
most common one in the iron pnictides, other mag¬ 
netic orders arise in Fe superconductors. FeTe orders 
in a double stripe pattern m and the alkaline iron se- 
lenides A^Fe 2 -a:Se 2 have block antiferromagnetism and 
a particular Fe vacancy arrangement [20l |2T]. Hole- 
doped Ba(Fei_a,Cra,)As 2 shows (7r,7r) order for x > 
0.3 [22]. Ferromagnetic order is much less usual than 
antiferromagnetism in Fe superconductors but it has 
been observed in a few cases UMl ESI M EMI 
1270] . [(Lii_a,Fea,)OH]FeSe based [T5UT^ and EuFe 2 As 2 
based compounds show FM (or canted AFM) 

in the hydroxide or Eu layer coexisting with superconduc¬ 
tivity in the EeSe or EeAs layer. EeTe has been observed 
to turn EM under pressure [SHE!]. This variety of mag¬ 
netic patterns indicates the magnetic order in these sys¬ 
tems is soft, namely, it can be affected by small changes 
in parameters. 

The magnetic softness has been addressed in the dif¬ 
ferent approaches used to study the magnetic state of 
these materials. In general, the magnetic ground state is 
dependent on the values of the parameters, chemical com¬ 
position and details of the structure m- Understanding 
the magnetic softness is important for superconductivity 
in these systems as (tt, 0 ) spin fluctuations lead to s~^ 
pairing symmetry |27l [80] while (tt, tt) fluctuations sup¬ 
press s~^ (and induce a d-wave state) |27Ij . EM correla¬ 
tions are expected to destroy singlet superconductivity. 


A. Magnetic softness in undoped compounds 

Weak coupling approaches assume the magnetic order 
is associated with a Eermi surface instability. Most of 
the Ee based superconductors show Eermi surface nesting 
between the F hole pockets and the X or Y electron 
pockets |3], see Fig. that would be expected to lead 
to (tt, 0) order. The alkaline Fe selenides, with no hole 
pocket in F and hence no Q = (tt, 0 ) nesting, do have a 
different (block AF) order [ 21 ]. However, some systems 
have a (tt, 0) nested Fermi surface but do not show a 
(tt, 0) magnetic ground state. This is the case of FeTe, 



U(eV) 


FIG. 8: Hartree-Fock phase diagram of a five-orbital 
model [ 31 ] as a function of the onsite interactions U and Jh for 
undoped (n = 6) Fe pnictides. The interorbital Coulomb re¬ 
pulsion U' is given by U' = U — 2Jh- Four different magnetic 
states have been considered. NM stands for non-magnetic, 
LM for low magnetic moment with violation of Hund’s rule 
(antiparallel orbital magnetic moments), and HM for a mag¬ 
netic state with parallel orbital magnetic moments. The (tt, tt) 
and double stripe regions (DS) also have a large value of the 
magnetic moment. Ferromagnetism (FM) is never the ground 
state for the parameters considered. Shaded areas correspond 
to insulating ground states (gap at the Fermi energy). Large 
JnjB promotes metallicity. The dashed lines and the labels 
on the right correspond to the results of calculating the en¬ 
ergies within a Heisenberg Hamiltonian for S — 2 spins with 
parameters calculated in second order perturbation theory for 
t/U < I. 


which orders in a double stripe pattern [TT] [with a Q = 
±(7r/2, ±7r/2)] which cannot be explained by nesting, and 
the phosphides, which do not order magnetically at all [TJ 

EH]. 

The magnetic ground state of a strong coupling model 
that includes the antiferromagnetic exchange to first Ji 
and second J2 neighbours is (tt, 0) for J1/J2 < 2 [33], and 
( 7 r, 7 r) for J 1 /J 2 > 2 . In order to get the double stripe 
pattern, the inclusion of exchange up to third neighbours 
J 3 is required because the number of antiferromagnetic 
and ferromagnetic links is equal for both first and second 
neighbours [ 86 ]. It has been argued that the stability 
of the double stripe order requires a very small Ji |272j 
(which is a likely scenario for relatively large values of 
Jh, see Fig. © or the addition of a biquadratic term 
K (2201 E23]. 

The relative success of the strong coupling approach to 
find the magnetic variability does not imply this model 
captures the full phenomenology of the Fe based super¬ 
conductors, as argued in Sec. in For instance, although 
the double stripe order of FeTe could be explained in 
terms of localised spins |22Qj , itinerant electrons were in¬ 
voked to explain the change in the effective spin per Fe 
from S' ^ I in the antiferromagnetic state to S ^ 3/2 
in the paramagnetic state revealed by an inelastic neu¬ 
tron scattering study of spin excitations |I93j . Similarly, 
the transition to a magnetic spiral with increasing y in 
Fei+^Te [TTJ 1274] could be explained by a modification 
of the exchange couplings through the RKKY interaction 
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between localised and delocalised electrons |275j . In gen¬ 
eral, models with localised and delocalised electrons have 
been proposed to explain the different magnetic phases as 
the competition between exchange couplings and kinetic 
energy driven order [71 13511361 [3^ [Tot]. 

Magnetic softness is also patent in ab-initio calcula¬ 
tions: different magnetic orders, with relatively close en¬ 
ergies, can be stabilised. Details in the chemical com¬ 
position and the structure may lead to different ground 
states, which are often associated to different magnitudes 
of the magnetic moment. The ground state for the pnic- 
tides is the (tt, 0 ) order, while for Fei+^Te it is a dou¬ 
ble stripe very closely followed in energy by a staggered 
dimer configuration [220] and (tt, 0) [ 86 ]. This proximity 
in energy is consistent with the coexistence of magnetic 
scattering at both Q = (tt, 0) and double stripe momenta 
reported in Fei.o 2 Tei_a,Sea, [276] . FeSe, which experi¬ 
mentally does not show a magnetic transition unless high 
pressures are applied [T 2 |, has a staggered dimer ground 
state within ab-initio calculations IMIEIZ]. Ferromag¬ 
netic order is usually not stabilised in DFT calculations 
for undoped systems [ 86 ], an exception being the 11 [ 220 ] 
and 111 compounds. Within DFT, in LiFeAs (which 
experimentally is a non-magnetic superconductor), both 
AFM and FM orders are possible, with AFM the ab-initio 
ground state. Calculations find a FM ground state in 
MgFeGe [278] and ferromagnetism in the related CuFeSb 
has been found in magnetisation and neutron scattering 
measurements [279] . These different magnetic orders - 
(7r,0), double stripe and FM - have been found to be 
the successive ground states for increasing values of the 
out-of-plane anion (the chalcogen or pnictogen) heights 
z. This relation has been explained from the different ef¬ 
fect of 2 : on the values of the exchange parameters, with 
J 3 getting stronger and Ji ^2 weaker as ^ increases, as 
calculated for FeSe^^Tei-a, [280] . 

The magnetic phase diagram of the undoped iron su¬ 
perconductors has been studied theoretically as a func¬ 
tion of the interactions (t/, and Jh) starting from dif¬ 
ferent five orbital model hamiltonians, as explained in 
Sec. [TT| Using a Hartree-Fock approximation in /c-space 
for the undoped (n = 6 ) iron-pnictides and comparing 
the energies of different magnetic orders [(7r,0), (7r,7r), 
double stripe, ferromagnetism, and non-magnetic (m = 
0)] for one of the models [31] results in the magnetic 
phase diagram m shown in Fig. Ferromagnetism is 
never found to be the ground state and the double stripe 
order that appears for large values of the interactions 
is accompanied by charge ordering (but charge order is 
not observed experimentally). The appearance of the 
double stripe order for large values of the interactions is 
consistent with the finding that chalcogenides are more 
strongly correlated than pnictides [98], see Sec. 

For intermediate values of the parameters, there is a 
close competition between the (7r,0) and (7r,7r) orders. 
This competition depends on the orbital filling (and it 
is hence related to orbital differentiation) and the crys¬ 
tal structure. A small variation of the crystal field can 
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FIG. 9: Strong coupling approach to the Fe-superconductors 
phase diagram. The exchange couplings Ji and J 2 are calcu¬ 
lated within second order perturbation theory from the tight- 
binding model in Ref. m assuming two different electron dis¬ 
tributions among the orbitals as sketched in (a). The x‘^ — 
configuration considers this orbital to be doubly occupied. 
Similarly for the configuration, (b) Exchange cou¬ 

plings to first Ji and second J 2 neighbours for the two pro¬ 
posed configurations. The vertical dotted lines separate the 
regions in which (7r,7r), (tt, 0) and FM are the ground states 
for increasing values of JhIU. (c) The exchange couplings 
calculated for half-filling (n = 5) always give rise to a (tt, tt) 
ground state. 


change the balance between these two orders: within 
Hartree-Fock, the transition between (tt, tt) and (tt, 0) 
is accompanied by a charge transfer from 32:^ — to 
[37]. Therefore, a slight shift down of the — y‘^ 
energy results in a more predominant (tt, 0) region on the 
phase diagram m This softness is a consequence of the 
very small crystal field splitting between the orbitals. In 
general, small changes in the orbital filling among differ¬ 
ent models affect the (tt, 7r)-(7r, 0) phase diagram [281] . 
For instance, the model in Ref. [89] has a relatively en¬ 
larged (7r,7r) region compared to the results shown in 
Fig.U 

The dependence of the ground state on the occupancy 
of the orbitals is also found in the strong coupling limit of 
the interacting models, see Eq. [^ Introducing the hop¬ 
ping as a second order perturbation to the interactions 
Hamiltonian, the five-orbital model can be mapped onto 
a Heisenberg model in which the couplings are calculated 
by adding the contributions from all the orbitals m- 
The exchange couplings to first Ji and second J 2 neigh¬ 
bours are calculated assuming that the orbitals are either 
singly or doubly (in the case of x‘^ — y^ and 3z‘^ — r^) 
occupied (namely, the magnetic moment corresponds to 
S = 2). Ji and J 2 can be of the same order due to the in¬ 
direct hoppings through the pnictogen or chalcogen m, 
which connect both first and second neighbours. 

Fig.i shows the resulting Ji and J 2 as a function 
of JhIU. Two orbital configurations are distinguished: 
the x‘^ — y^ configuration and the 3z‘^ — configuration, 
whose name refers to which orbital is doubly occupied, 
see Fig. |^(a). The values of the exchange couplings Ji ^2 
depend on the orbital configuration, the hoppings, and 
the interactions. The vertical lines in Fig.[^(b) separate 
the regions with different magnetic order. These regions 
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are superimposed with dashed lines in the Hartree-Fock 
phase diagram in Fig. An increasing Hund’s coupling 
Jh suppresses both Ji and J 2 [371142]. Ji decreases much 
faster than J2, such that the condition for (tt, 0) order 
(Ji < 2J2) is fulfilled for a certain range of Jh values. 
Consistent with Hartree-Fock results, the (tt, 0) order cor¬ 
responds to the —y^ configuration and the (tt, tt) order 
to the 3z^ — configuration. This relation arises from 
the relatively large value of the x‘^ — y‘^ hoppings between 
first neighbours m and reveals the delicate balance be¬ 
tween Jh and the crystal field. The largest values of Jh 
change the sign of the exchange couplings leading to FM 
order. The suppression of Ji and J 2 is favoured by vir¬ 
tual transitions involving the doubly occupied 3z^ — 
or x‘^ — y^ orbitals [37] . 

Ji and J 2 also depend on the angle a that charac¬ 
terises the Fe-As tetrahedra, see Fig.[^ through the value 
of the hoppings m leading to a predominance of (tt, 0) 
for squeezed tetrahedra and FM for elongated tetrahe¬ 
dra [^ . 



FIG. 10: Hartree-Fock phase diagram as a function of doping 
with the same multi-orbital model as considered in Fig. [^ 
The possible magnetic orders considered are (7r,0), (7r,7r), 
double stripe and FM. The value of the Hund’s coupling is 
fixed to Jh = 0.22U. Hole doping (n < 6) promotes anti¬ 
ferromagnetic interaction and the (tt, tt) order is found close 
to n = 5. Electron doping n > 6 promotes the double stripe 
phase. There are big regions of phase separation in the phase 
diagram. 


B. Magnetic softness in doped compounds 

The magnetic softness also appears upon doping. At 
the extreme case of n = 5 [Fig.|^(c)], the strong coupling 
limit gives Ji / J2 > 2 for all values of Jh and the system 
has (tt, tt) order. The competition between different mag¬ 
netic orders as a function of the filling n has been studied 
within the Hartree-Fock approximation in the recipro¬ 
cal and in the real |lQ8j spaces. The general trends 
of both sets of results are that (tt, tt) order is the ground 
state around n = 5, which is expected at half-filling where 
correlations are stronger and the strong coupling result 
applies, and that the tendency towards FM increases with 
electron doping, with a (tt, 0) region around n = 6, see 
Fig.[^ Phase separation appears in different parts of the 
phase diagrams. Many-variable variational Monte Carlo 
calculations for ab-initio low-energy models find the same 
qualitative features as a function of doping |114j . 

These theoretical results are consistent with measure¬ 
ments on materials in which Fe is substituted by a neigh¬ 
bouring transition metal while keeping the same crystal 
structure. For instance, Mn^+ ions are half-filled (d^) 
while Co^+ has 7 electrons (dJ). In BaMn2As2 |282j 
or BaMnPnF (Pn=As, Sb, Bi) [283], a (7r,7r) AF, in¬ 
sulating phase is found with very high T/v and magnetic 
moment. On the electron-doped side of the phase di¬ 
agram, LaCo(As,P)0 are itinerant ferromagnets [284], 
CaCo2As2 has FM planes [285], while BaCo2As2 is para¬ 
magnetic but allegedly close to ferromagnetism |286j . 
From this, we should expect a tendency towards (7r,7r) 
order when doping with holes and towards ferromag¬ 
netism when doping with electrons. Consistently, neu¬ 
tron diffraction experiments [22] in Cr (d^) doped sam¬ 
ples Ba(Fei_a,Cra,)As2 show a transition between the 
(tt, 0) and the (tt, tt) order for x ^ 0.3. Note that partially 
substituting Fe by Cr or Mn does not produce a transi¬ 


tion to a superconducting phase |287j . It is believed that 
Cr and Mn do not incorporate itinerant electrons in the 
samples but act as magnetic impurities. 

Neutron and X-ray diffraction measurements on 
Ba(Fei_a,Mna,)2As2 reported the suppression of the or¬ 
thorhombic to tetragonal structural transition at x ^ 0.1 
but that the Bragg peaks corresponding to stripe order 
[(tt, 0) or (0, tt)] survive in the tetragonal lattice for x > 
0.1 |288j . Inelastic neutron scattering experiments |289j 
reported the coexistence of both (tt, 0) and (tt, tt) spin 
fluctuations for Ba(Feo.925Mno.075)2As2. Fluctuations at 
(tt, 0) have a strong temperature dependence with its sig¬ 
nal becoming broader above T^v- In contrast, (7r,7r) fluc¬ 
tuations are very short ranged and weakly T dependent, 
they are not suppressed when the (tt, 0) long range order 
sets in, and persist up to room temperature [289]. T^v is 
suppressed upon Mn substitution up to x ^ 0.1 and it is 
enhanced for larger x while the transition is significantly 
broadened |287[ I288j . However, neutron Larmor diffrac¬ 
tion together with nuclear magnetic resonance, muon 
spin-relaxation and neutron diffraction experiments |290j 
for X > 0.1 samples point to the coexistence of (tt, 0) or¬ 
dered orthorhombic areas and paramagnetic tetragonal 
areas. Antiferromagnetic short range ordered regions are 
seen to form at temperatures much higher than T/v |29Qj . 
This interpretation as a coexistence phase is consistent 
with a model calculation that includes the interaction 
between the conduction electrons and the Mn magnetic 
impurities |291j . 

A tetragonal (C4 symmetry) phase has also been seen 
in various 122 systems close to the suppression of the 
antiferromagnetic phase, where there is a magnetism- 
superconductivity coexistence region. So far three com¬ 
pounds have shown this reentrant tetragonal phase: 
Bai_a,Naa,Fe2As2 |292j . Bai_a,Ka,Fe2As2 |293j . and 
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Sri_a^Naa:Fe 2 As 2 |294j . Bai_a^Naa^Fe 2 As 2 shows a C 4 
phase coexistent with the orthorhombic one just at 
the antiferromagnetic transition suppression for a small 
range of doping |292j , while this tetragonal region is fully 
inside the C 2 one in Sri_a^Naa^Fe 2 As 2 |294j . The tetrag¬ 
onal phase in Bai_a^Ka;Fe 2 As 2 is also at the magnetic 
transition but, differently to the other two cases, it is 
also reentrant as a function of temperature, being sup¬ 
pressed as the superconductivity sets in |293j . The reen¬ 
trant phase is accompanied by a spin reorientation |292j . 
Polarised neutron diffraction data in single crystal sam¬ 
ples of Bai_a,Naa,Fe 2 As 2 have identified the spin reorien¬ 
tation as occurring towards the out-of-plane c-axis |295j . 

The magnetic order in the reentrant tetragonal C 4 
phase is a matter of active research at the moment. One 
possibility, consistent with a nesting origin of magnetism 
is the so-called double-Q order |123[ 1245112881129211296F 
I298j . The double-Q order is a combination of ( 7 r, 0 ) and 
( 0 , 7 r) which fulfils the C 4 symmetry. This phase has 
been found within weak coupling models |123j . on a two 
orbital |245j and on a five orbital interacting model with¬ 
out disorder |298j within unrestricted Hartree-Fock, and, 
for the case of Mn-doped BaFe 2 As 2 , it has been argued 
via Ginzburg-Landau calculations to arise as the ground 
state between the (tt, 0) C 2 phase and the (tt, tt) order 
expected when the local moments on Mn dominate |296j . 
However, different magnetic and orbital order configu¬ 
rations may be compatible with the C 4 structure. It 
has been argued that a full characterisation of this phase 
could gives clues on the determination of the magnetic or 
orbital origin of the magnetoestructural transition |299j . 


VI. SUMMARY AND CONCLUSIONS 

The magnetic phase diagram of iron superconductors 
gives us information about the underlying electronic cor¬ 
relations in these systems. Such information may be cru¬ 
cial to understand the origin of the high Tc supercon¬ 
ductivity. Here we have reviewed different proposals to 
explain the magnetic interactions in iron superconduc¬ 
tors. The multiorbital character of these systems under¬ 
lies their complex phenomenology and must be consid¬ 
ered for the correct description of the electronic proper¬ 
ties and the comparison with experiments. The magnetic 
interactions are directly linked to the strength of corre¬ 
lations in these materials, with Hand’s coupling playing 
a very important role. Taking into account the orbital 
degree of freedom is crucial (a) to unify, in a single U-Jh 
phase diagram, the main mechanisms proposed to explain 
the (tt, 0 ) magnetic order; (b) to explain the magnetic 
softness; (c) to address doping dependent properties; and 
(d) to understand the anisotropies. 

Three different regions are distinguished in the quasi¬ 
particle weight phase diagram of iron pnictides plotted in 
Fig. (i) a metallic region with > 0.5 for moderate 
interactions, (ii) a strongly correlated Hund metal with 
small quasiparticle weight, spin-polarized local atomic 


states, and orbital differentiation (i.e., orbital dependent 
correlations), for intermediate U and Jh^ and (iii) a Mott 
insulating state at large U. The metallic regions (i) and 
(ii) are separated by a crossover. The interaction b^Hund 
at which this crossover takes place decreases in hole- 
doped materials, as the occupation of the Fe-d orbitals 
approaches half-filling. 

We argue that all the Fe-superconductors can be placed 
within this unified phase diagram. While the strength of 
correlations is still under debate, theoretical estimates 
and comparison with experiments suggest that the iron 
arsenides have JrIU ^0.15 — 0.3, are not far from the 
crossover between the two metallic states, and become 
more correlated with hole-doping. Iron chalcogenides 
are believed to be more correlated, but still in the Hund 
metal region, or, in the case of the alkaline-doped chalco¬ 
genides, close to the Mott insulating state. 

The different states which appear in the (tt, 0) Hartree- 
Fock phase diagram (Fig.|^ can be related to those in the 
correlations phase diagram (Fig. |^, and to weakly and 
strongly coupling descriptions proposed within different 
approaches, see discussion in Sec. |mg While Hartree- 
Fock cannot address the localisation of quasiparticles, the 
behavior and properties of the insulating state at large 
U compares well with expectations from the strong cou¬ 
pling Heisenberg approach and the Mott insulating state. 
The magnetic region with strong orbital differentiation in 
Fig. seems related to the Hund metal state in Fig. 

In this magnetic region zx^ 3^^ — and x‘^ — are itin¬ 
erant, while yz and xy are half-filled and gapped. The 
interplay between these two sets of orbitals suggests a 
double-exchange mechanism for magnetism with itiner¬ 
ant and localized orbitals. Finally, the metallic region 
with moderate correlations in Fig. is connected with 
three regions in the (tt, 0 ) magnetic phase diagram: a 
non-magnetic state at small U, a magnetic ordered state 
which violates Hand’s coupling, and a nesting originated 
itinerant magnetic state. Within this (tt, 0) magnetic 
phase diagram, the Fe arsenides are most probably within 
the crossover between the itinerant and the orbital dif¬ 
ferentiated regions. 

The multiorbital character of iron superconductors al¬ 
lows for different magnetic ground states to arise in un¬ 
doped systems. For intermediate values of U and Jh^ 
( 7 r, 0 ) and ( 7 r, 7 r) compete. This competition can be 
tuned with a slight change in the relative crystal field 
between x‘^ — y^ and 3z‘^ — r^. Larger values of the in¬ 
teractions, specially Hand’s coupling, stabilise the dou¬ 
ble stripe, which is precisely the magnetic order found 
in iron chalcogenides, with strong correlations. Hole 
doping enhances correlations taking the systems closer 
to a Mott insulator and a ( 7 r, 7 r) order, while electron¬ 
doping promotes itinerancy and a tendency towards fer¬ 
romagnetism, in accordance with experimental observa¬ 
tion. Different magnetic states can be stabilised through 
small changes in the crystal structure, like the pnictogen 
(or chalcogen) height. 

The C 4 symmetry breaking of the (tt, 0) antiferromag- 
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netic state splits the degeneracy of the yz and zx or¬ 
bitals and gives rise to orbital ordering as corroborated 
in multiorbital models. Spontaneous generation of or¬ 
bital ordering has not been found without magnetism in 
realistic five orbital models. These models also find that 
both orbital ordering and magnetic moment are anticor¬ 
related with the Drude anisotropy in the magnetic state. 
The topology and the morphology of the Fermi surface 
affects strongly to the resistivity anisotropy in the mag¬ 
netic state in presence and in absence of impurities. The 
tetragonal symmetry breaking of magnetism is also be¬ 
hind the anomalous Raman response of the Aig phonon 
when undergoing the magnetic transition. 

Effective models accounting for the orbital degree of 
freedom are necessary to understand the nematic phase. 
Spin fluctuations give rise to dxzj dyz orbital splitting 
in this phase. Spin and orbital fluctuations are present 
then in the nematic state even in the spin-nematic sce¬ 


nario. The anisotropic properties of the spin fluctuations 
in the nematic phase depend sensitively on the singular 
C 4 symmetry of the yz and zx orbitals. 

Seven years of intense research on iron superconduc¬ 
tors have given rise to an unparalleled development of 
new concepts involving correlations in multi-orbital sys¬ 
tems. Understanding their implications is likely the key 
to elucidate the mechanism of high Tc superconductiv¬ 
ity in these compounds. Iron superconductors promise 
to keep theorists and experimentalists busy for years to 
come. 
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